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ABSTRACT 


This document Is a description of and users manual for 
a USA FORTRAN IV computer program which determines 
boundary conditions for a thin wing lifting surface 
program. This program, the geometry program, and 
several other programs are used together in the 
analysis of lifting, thin wings in steady, subsonic 
flow according to a kernel function lifting surface 
theory. The program calculates specific types' of 
boundary conditions completely automatically such as 
those necessary to determine pitch and roll damping 
derivatives. The program also accepts descriptions of 
the camber or downwash and twist in the form of tables 
and/or coefficients of equations. The program performs 
Interpolations so that tables and/or coefficients can 
apply at stations selected by the user and not at 
stations dictated by the control point locations. The 
program uses Information stored on a geometry file and, 
optionally, on an influence matrix file. The boundary 
conditions that it calculates are stored on a boundary 
condition file. An equation solving program will read 
the influence matrix and boundary condition files and 
determine the coefficients in the expansion for the 
lifting pressure distribution. 
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1 INTRODUCTION 


This document is a description of and users manual for 
a USA FORTRAN IV computer program which . determines 
boundary conditions for a thin wing lifting surface 
program. This program, the geometry program, and 
several other programs are used together in the 
analysis of lifting, thin wings in steady, subsonic 
flow according to a kernel function lifting surface 
theory. The program calculates specific types of 

boundary conditions completely automatically such as 
those necessary to determine pitch and roll damping 
derivatives. The program also accepts descriptions of 
the camber or downwash and twist In the form of tables 
and/or coefficients of equations. The program performs 
interpolations so that tables and/or coefficients can 
apply at stations selected by the user and not at 
stations dictated by the control point locations. The 
program uses Information stored on a geometry file and, 
optionally, on an influence matrix file. The boundary 
conditions that It calculates are stored on a boundary 
condition file. An equation solving program will read 
the influence matrix and boundary condition files and 
determine the coefficients In the expansion for the 
lifting pressure distribution. 

Questions concerning either this document or the 
compute r program or the associated computer programs 
should be directed to: 

R. T. Medan 

Mall Stop 247-1 

Ames Research Center 

Moffett Field, Ca., 94035 
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2 USER'S INSTRUCTIONS 

2.1 INITIAL SETUP FOR AMES' TSS SYSTEM 

For either batch or conversational processing the 
following TSS commands must be given. These commands 
are required once and only once for each user ID. 
The first three commands create the identification 
number file named IDFILE. This file contains four 
zeroes in binary form. 

SHARE MEDAN, FSARTM, INIHFILE 
CDS MEDAN, IDFILE 
DELETE MEDAN 

SHARE MEDAN, FSARTM, LSPROG. VI 

2.2 CONVERSATIONAL USE ON AMES' TSS SYSTEM 

All integer data should be entered in a 16 f 5 
format, all floating point data in 8F10.0 
format, and all logical data in 10L1 format. 
See section 8 for an example of a terminal 
sess 1 on . 

USER: After logging on enter the following: 

AMES USYSLIB 
JOBLIBS SYSULIB 
JBLB MEDAN 

It isn’t necessary to issue DDEFs for anything 
except the input data since the program Issues 
them using the subroutines AIMFIL, GEMFII, and 
BCFIL. 

USER: CALL BC$ 

PROG: ENTER BATCH 

USER: Enter carriage return for conversational mode. 

PROG: ENTER INPUT DEVICE NUMRER 

USER: Hit carriage return for terminal input. 

Otherwise input starting with (PCS) will be 
from a dataset on disk as referenced by the 
input device number. This number must be a 
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positive Integer from 1 through 99 excluding 
4, 5 , 6, 7, 8, 9, and XI. This dataset must 

be named In a DDEF statement or else AMES' TSS 
will expect terminal Input. If this method of 
Input Is used, then the program will not give 
conversational prompts for the data starting 
with (BCS). 

To terminate execution, enter a negative 
numbe r. 

PROG: ENTER ODISK 

USER; For terminal output enter carriage return. 

For output to a disk file enter a positive 
non-zero number less than 10. For TSS the 
output will be found on the file named 
OUTPUT. BC .NJJ, where X is the numerical value of 
ODISK. The program Issues its own DDEF 
conmands for the output file so no control 
cards are needed. This output Is written on 
logical un 1 t 4 . 

PROG: ENTER ID1,ID2 

USER: Enter Identification numbers. 

I D1-t I dent I f ica 1 1 on number of the geometry 
file. If a negative number is entered, then 
the most recent geometry file will be used. 
Enter zero to terminate execution. 

ID2 — Identification number of AIM file. Enter 
a nonzero number only If the control points 

have been changed by the influence matrix 

program. If a negative number is entered, 
then the most recent AIM file will be used. 
Enter zero if no AIM file has been created or 
if the control points have not been changed by 
the influence matrix program. 

PROG: ENTER { BCS)/{ BCAS)/PR I NTA, PR I NTB,NSTORE 

USER: Enter the array (BCS), carriage return. 

Enter the array (BCAS), carriage return. 

Enter PRINTA, PRINTB, NSTORE . 
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Each element of (BCS) and (BCAS) corresponds 
to a different downwash condition. A true 
value means that the corresponding boundary 
condition will be calculated and, if NSTORE is 
.FALSE., then written on a file for use in the 
equation solving program. Both arrays are 
read and acted upon tf the wing is symmetric 
or unsymmetric. In the latter case the 
boundary condition will be calculated for the 
entire wing. See section 3 for a description 
of the types of cases referenced by (BCS) and 
(BCAS) and for auxilliary inout required for 
each element of (BCS) and (BCAS). 

PRINTA is a logical variable which should be 
set to .TRUE, if all the downwash modes which 
are calculated should also be printed 
(regardless of the value of PRINTB). PRIMTB 
is a logical variable which should be set to 
.TRUE. if the 5th-10th symmetric and the 
3rd-10th unsymmetric modes which are 
calculated should also be printed. NSTORF is 
a logical variable which should be set to 
.TRUE, if the boundary conditions should be 
calculated but not stored. 

PROG: ENTER LL (if there is a flap or flaps) 

USER: Enter LL. 

LL is the number of times that subroutine 
FLPDWN will intergrate to find the flap 
downwash mode or modes. Successive 
integrations will use more points than 
previous ones up to a maximum of JJMAX, the 
maximum number of integration stations 
available from the geometry file. Only the 
results for the largest number of Integration 
points will be retained and written on the 
boundary condition file. 

PROG: The program will ask for auxilliary 

information (if any) required for the various 
cases referred to by (BCS) and (BCAS). See 
section 3 for the auxilliary information 
required. After this is done, the program 
will loop back to the point in the main 
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program at which the Input device number Is 
requested. The user can run another case at 
this point or else enter a negative number to 
terminate execution. 

2.3 ' AMES’ TSS BATCH JOBS 

The batch mode operates the same as the 

conversational mode with the following exceptions: 

(1) A "T 11 should be put In column 1 of the first 
card. 

(2) The Input device number must not be entered 
since all Input in this case Is assumed to originate 
on unit 5. 

(3) ODISK must not be entered since the program 
assumes all output should go onto unit 6. 

(4) To terminate execution In the batch mode, enter 
zero for IP1 (i.e. a blank card at the end of the 
input deck of the final case). 

2.4 CONVERSION TO OTHER COMPUTERS 

Remove all calls to GEMFIL, BCFIL, AIMFI1, OREY, and 
CVRT and use appropriate control cards. These 
subroutines Issue DDEFs and RELEASE commands making 
control cards unnecessary on TSS. Only the main 
program needs to be changed. These, hopefully, are 
the only changes that need to be made since 

considerable effort was made to program In standard 
FORTRAN. 
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3 DOWNWASH DISTRIBUTIONS 

In the following explanations of the downwash modes, 
XSI and ETA are chordwise and streamwise coordinates, 
respectively, normalized by the effective semi-span 
(the semispan after yawing). Any Input required is 
prompted for in the conversational mode, if input is 
from the terminal. 

3.1 BCS(l) uniform mode 
ALFA(XS I , ETA) “1 . 

INPUT REQUIRED: NONE 

3.2 BCS( 2) pitching mode 

ALFACXS I , ETA) =(XS l/BRAT-XS I CM)/CBARRR where 

XS I CM=center of mass or other reference position. 
XSICM is to be given in the coordinate system fixed 
to the wing while XSI is in the wind-centered 
coordinate system This mode corresponds to a 
non-dimensional pitch rate of 1. The solution for 
this mode gives the wing contribution to the Q 
stability derivatives (in the wind centered 
coordinate system). BRAT is the ratio of effective 
to actual semi span (ref. 1). 

INPUT REQUIRED: XSI CM 

3.3 BCS(3) linear, symmetric twist mode 
ALFA(XSI , ETA) =ABS( ETA) /BRAT 

INPUT REQUIRED: NONE 

3.4 BCS(4) parabolic, symmetric twist mode 
ALFA( XSI, ETA) =( ETA/BRAT) **2 

INPUT REQUIRED: NONE 

3.5 BCS(5) residual flap downwash mode. 

If there are 2 flaps, then this mode Is for the 
symmetric deflection of both. Refer to NASA TN 
D-7251 and to subroutine FLPDWN for further 
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documentation. The flap downwash modes are the only 
ones affected by the Mach number. The value of 

BCS(5 ) and BCAS(3) will be ignored if there is no 
flap data available from the geometry file, 

3.6 BCS(6)-BCS(10) 

These are modes for which tables and/or certain types 
of coefficients will be read by the program. The 
modes are defined by streamwise distributions 
(coefficients or tables) at spanwise stations 
selected by the user. These modes are assumed 
symmetric only If the wing is symmetric. See 
subroutine BOUND for further documentation. 

INPUT REQUIRED: SEE SECTION 4 

3.7 BCAS(l) anti-symmetric uniform mode (zero at 
center). 

INPUT REQUIRED: NONE 

3.8 BCAS(2) rolling mode (right tip down) 

ALFA( XS I , ETA) ®ETA/ BRAT 

This mode corresponds to a non-d imens i onal rol 1 rate 
of 1.0. The solution for this mode gives the wing 
contribution to the quasi-static, P stability 

derivatives (in the wind centered coordinate system). 

INPUT REQUIRED: NONE 

3.9 BCAS(3) anti-symmetrlc residual flap mode. 

If the wing is unsymmetric, there is no need to 
invoke this mode. 

3.10 BCAS( 4 ) - BCAS ( 10 ) 

For a symmetric planform these modes are the 

anti -symmetric counterpart of those referred to by 
BCS(6)-BCS(10) . For an unsymmetric wing the 

treatment of these modes is the same as those 
corresponding to BCS(6)-RCS(10) . 

INPUT REQUIRED: SEE SECTION 4 
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4 INPUT FOR SUBROUTINE SLOPES 

All Integer data should be entered in a 1615 format 
and all floating point data in 8F10.0 format. 

4.1 The following Input items are requested each time 
SUBROUTINE SLOPES Is called/ i.e. for each true value 
of BCS(6)-BCS( 10) and BCAS(4)-BCAS(10) . See section 
8 for examples on the use of subroutine SLOPES. 

PROG: ENTER NSPSEC, ITYPES 

USER: Enter NSPSEC, (TYPES. 

NSPSEC is the number of spanwise sections at 
which twist and/or camber data will be 
supplied by the user. 

(TYPES is the type of spanwise interpolation. 

0 Implies straight line interpolation. 

+0 implies that C0D1M2 will be used(ref. 1). 

PROG: ENTER XK (only if ITYPES+0). 

USER: Enter XK. 

This Is an Interpolation control constant used 
for the spanwise interpolation when ITYPES+0. 
A value of 0. will give linear interpolation 
in the end intervals. A value of 1.0 will 
give parabolic Interpolation in the end 
intervals. A value in between will give a 
curve in between. 

4.2 The data in this section will be requested a total of 
NSPSEC times In the order given (i.e. ETA/ TWIST; 
ITYPEC; data required for specific value of ITYPEC; 
ETA/ TWIST; etc. repeated a total of NSPSEC times). 

PROG: ENTER ETA/ TWIST 

USER: Enter ETA/ TWIST. 

ETA is the spanwise station at which the twist 
and camber data apply. It may exceed 1 in 
value in order to control the C0DIM2 
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Interpolation near the ends If C0DIM2 Is used. 

TWIST Is the angle of twist In radians. 
Positive twist will tend to increase the angle 
of attack of the wing. 

PROG: ENTER ITYPEC 

USER: Enter ITYPEC. 

ITYPEC denotes the type of chordwlse slope 
definition that the user will use. 

4.2.1 t TYPEC=0 

In this case there is no camber, just twist. 
INPUT REQUIRED: None. 

4.2.2 I TYPEC=1 

In this case a set of polynomials In CHI 
defines the slope distribution, whtch is the 
derivative of z/chord with respect to CHI. 
CHI is the local chordwise variable such that 
0<CHI<1. 

INPUT REQUIRED: NPOLS/CH I MAX, SCALE, CH I REF, 

( NPOLS t l mes)/ POLYNOMIAL COEFFICIENTS, 

( NPOLS times). 

NPOLS is the number of polynomials. CHIMAX 
denotes the upper limit of the polynomial. 
The default value is 1. SCALE is a scale 
factor for the polynomial coefficients. The 
default value is 1. CHIREF Is the origin of 
the polynomial. The degree of each polynomial 
Is limited to 8 or less. This data is 
prompted for and entered in the following 
manner: 

PROG: ENTER NPOLS 

USER: Enter NPOLS. 

PROG: ENTER CHI MAX, SCALE, CH I REF/ 

POLYNOMIAL COEFFICIENTS 
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USER; Enter CHIMAX# SCALE# CHIREF# carriage return. 
Enter POLYNOMIAL COEFFICIENTS. 

CHIMAX# SCALE# CHIREF# and the polynomial 
coefficients will be requested a total of 
NPOLS times. 

4.2.3 I TYPE C= 2 

In this case a set of polynomials in CHI 
defines the camber distribution (z/chord as a 
function of CHI). This set of polynomials 
will be differentiated to determine the slope, 
so the leading coefficient in each polynomial 
is actually Irrelevant. 

INPUT REQUIRED; NPOLS/CH I MAX# SCALE# CH I RFF# 
(NPOLS t imes)/ POLYNOMIAL COEFFICIENTS# 

(NPOLS times). 

NPOLS Is the number of polynomials. CHI MAX 
denotes the upper limit of the polynomial. 
The default value is 1. SCALE Is a scale 
factor for the polynomial coefficients. The 
default value is 1. CHIREF is the origin of 
the polynomial. The degree of each polynomial 
is limited to 8 or less. This data is 
prompted for and entered in the following 
manner: 

PROG: ENTER NPOLS 

USER: Enter NPOLS. 

PROG; ENTER CH I MAX, SCALE# CH I REF/ 

POLYNOMIAL COEFFICIENTS 

USER: Enter CHIMAX# SCALE# CHIREF# carriage return. 

Enter POLYNOMIAL COEFFICIENTS. 

CHIMAX# SCALE# CHIREF# and the polynomial 
coefficients will be requested a total of 
NPOLS times. 

4.2.4 ITYPEC=3 
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In this case a table of values defines the 
slope dlstrlbut Ion, the derivative of z/chord 
with respect to CHI. Subroutine C0DIM2 will 
be used for the interpolation. 

INPUT REQUIRED: SCALE, XKC/ table of (CHI, 

ALFA) pa i rs . 

SCALE is a scale factor for the table entries. 
The default value Is 1. XKC Is the end point 
interpolation control for C0DIM2. If XKC=0., 
then straight lines will be used in the end 
intervals. If XKC=1., then full parabolic 

interpolation will be used in the end 
intervals. A value In between will give a 
curve in between. CHI Is the local chordwise 
variable such that 0<CHI<1 (values outside of 
this range may, however, be entered in the 
table to control the interpolation near the 
end points). ALFA is the derivative of 

z/chord with respect to CHI. A maximum of 40 
.table entries are allowed not including 

CH 12.99.0, which Is the value used to mark the 
end of the table. This data is prompted for 
and entered in the following manner: 

PROG: ENTER TABLE SCALE FACTOR, A1(0R XKC), 

AN( OR EPS), AND THE TABLE, ( CHI =99 .STOPS) 

USER: Enter SCALE, XKC, carriage return. 

Enter the (CHI, ALFA) table. 

CHI and ALFA are pairs In the table. Each 
pair goes on a separate line or card. 

CHI2.99.0 marks the end of the table. 

4.2.5 ITYPEC=4 

In this case a table of values defines the 
camber distribution. C0DIM2 will be used to 
determine the camber in the vicinity of the 
control points. Then numerical 

differentiation will be used to determine the 
slopes . 

INPUT REQUIRED: SCALE, XKC, EPS/ table of 

(CHI, z/chord) values. 
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SCALE is a scale factor for the -table entries* 
The default value is 1. XKC is the end point 
interpolation control for C0DIM2. If XKC=0., 
then straight lines will be used In the end 
intervals. If XKC=1., then full parabolic 
interpolation will be used in the end 
intervals. A value in between will give a 
curve in between. EPS is used for numerical 
differentiation of the camber distribution 
determined by a table and C0PIM2 (controlled 
deviation interpolation method). The default 
value is .005. CHI is the local chordwise 
variable such that 0<CHI<.1 (values outside of 
this range may, however, be entered in the 
table to control the Interpolation near the 
end points). A maximum of L0 table entries 
are allowed not including CHI >9 9.0, which is 
the value used to mark the end of the table. 
This data is prompted for and entered In the 
fol lowing manner: 

PROG: ENTER TABLE SCALE FACTOR, Al( OR XKC), 

AN( OR EPS), AND THE TABLE, (CH 1=99, STOPS) 

USER: Enter SCALE, XKC, EPS, carriage return. 

Enter the (CHI, z/c) table. 

CHI and z/c are pairs in the table. Each pair 
goes on a separate line or card. CHt>99.0 
marks the end of the table. 

4.2.6 ITYPEC=5 

In this case a table of values defines the 
slope distribution. Cubic spline fits will be 
used to determine values at the chordwise 
control stations. Two types of spline fits 
are used. In the first (subroutine SPLINl) 
the end point derivatives are not used. In 
the other (SPLIN2) the end point derivatives 
are used. The latter program generally gives 
a better fit. 

INPUT REQUIRED: SCALE, Al, A2/table of (CHI, 

ALFA) values. 
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SCALE Is a scale factor for the table entries. 
The default value Is 1. A1 and AN are the 

derivatives of the slope at the end points of 
the table. A1 and AN should be put In If 
possible. If both values are actually zero, 
then input one as l.E-30. If either Is given, 
then both must be given. If both A1 and A2 
are zero then SPLINl will be used. Otherwise 
SPLIN2 will be used. CHI Is the local 

chordwise variable such that 0<CH!<1 (values 
outside of this range may, however, be entered 
In the table to control the Interpolation near 
the end points). ALFA is the derivative of 
z/chord with respect to CHI. A maximum of 40 
table entries are allowed not including 
CHI>99.0, which Is the value used to mark the 
end of the table. This data Is prompted for 
and entered In the following manner: 

PROG: ENTER TABLE SCALE FACT0R,A1(0R XKC), 

AN(OR EPS), ANP THE TABL E( CH I *99 . STOPS ) 

USER: Enter SCALE, Al, AN, carriage return. 

Enter table of (CHI, ALFA) values. 

CHI and ALFA are pairs In the table. Each 
pair goes on a separate line or card. 
CHI>99.0 marks the end of the table. 

4.2.7 ITYPEC=6 

In this case a table of values defines the 
camber distribution. A cubic spline fit will 
be determined for the table and then 
differentiated to obtain the slopes. Two 
types of spline fits are used. In the first 
(subroutine SPLINl) the end point derivatives 
are not used. In the other (SPLIN2) the end 
point derivatives are used. The latter 
program generally gives a better fit. 

INPUT REQUIRED: SCALE, Al, A2/tab!e of (CHI, 

z/chord) values. 

SCALE Is a scale factor for the table entries. 
The default value is 1. Al and AN are the 
slopes at the end points of the tahle. Al and 


13 



Boundary Condition Program 


AN should be put In if possible. If both 
values are actually zero, then input one as 
l.E-30. If either is given, then both must be 
given. If both A1 and A2 are zero then SPLIM1 
will be used. Otherwise SPLIN2 will be used. 
CHI is the local chordwise variable such that 
(KCHI.<1 (values outside of this range may, 
however, be entered in the table to control 
the interpolation near the end points). A 
maximum of 40 table entries are allowed not 
including CHI>99.0, which is the value used to 
mark the end of the table. This data is 
prompted for and entered in the following 
manner: 

PROG: ENTER TABLE SCALE FACTOR, A1C0R XKC), 

AN( OR EPS), AND THE TABLE, ( CH I =99 .STOPS ) 

USER: Enter SCALE, Al, AN, carriage return. 

Enter table of (CHI, ALFA) values. 

CHI and ALFA are pairs in the table. Each 
pair goes on a separate line or card. 
CHI>99.0 marks the end of the table. 

4.2.8 ITYPEC>6 or ITYPEC<0 

This causes SUBROUTINE USLOPE to be called, 
which is a program to be supplied by the user 
in the situation that he wants to define the 
chordwise slope distribution at some 
particular spanwise station by a method 
unavailable in subroutine SLOPES. See section 
5.4 for an explanation of the arguments to 
subroutine USLOPE. 
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5 PROGRAM DESCRIPTIONS 

5.1 MAIN PROGRAM 

This program determines boundary conditions for a 
thin wing lifting surface program. It uses 
Information stored on the geometry file and, 

optionally, on the influence matrix file. The 

boundary conditions that It calculates are stored on 
a boundary condition file. An equation solving 
program will read the influence matrix and boundary 
condition files and determine the coefficients in the 
expansion for delta-c.p.. 

This program reads from the geometry file (unit 7), 
the influence matrix file (unit 11), and the 
identification number file (unit 20). It writes on 
the identification number file and the boundary 
condition file (unit 8). See section 2 and 3 for 
conversational and batch use. 

Arrays ETA, STHETA, TANLEL, TANLER, TANTEL, TANTFR, 
CORDIP, XSILIP, CHIFPI, C2IP should be dimensioned at 
least as large as JJMAX. Arrays N INDEX and ETACP 
should be dimensioned at least as large as MM. ALFA 
should be dimensioned as large as the number of 
control points being used. CHICP should be 
dimensioned as large as PP. 

5.2 SUBROUTINE BOUND 

This subroutine is called by the main program to 
determine downwash cases. 

This subroutine writes on the boundary condition file 
(unit 8), See section 2 for conversational and batch 
use . 

INPUT VARIABLES: 

PP Number of CHICP values (i.e. the number of 

chordwise control points) 

(CHICP) An array containing chordwise stations at 
which boundary conditions will be 
determined. 

MMP Number of spanwise control points on wing 

for symmetric boundary conditions. 
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MM PA 

( ETACP) 

CBARBR 

BRAT 
( BCS) 

( BCAS) 

FLAP 

PRINTA 

PRINTB 

NSTORE 

SYMM 

CONV 

C0NV2 

UCI 

UCO 

U5 

U6 

U8 

JRATIO 

LL 


Number of spanwlse control points on wing 
for antisymmetric boundary conditions. 
Spanwlse stations at which boundary 
conditions will be determined. 

Ratio of chordwlse reference length to 
effective semispan. 

The ratio of effective to actual semispan. 

An array of downwash condition Indicators. 
See section 3. 

An array of downwash condition indicators. 
See section 3. 

Logical variable, whose value Is to be 
.TRUE, if there is a control surface (the 
program currently does not support wings 
with control surfaces, however). 

Logical variable whose value is .TRUE, if 
all the downwash modes which are calculated 
should be printed (regardless of the value 
of PRINTB). 

Logical variable whose value is .TRUE, if 
the 5th-10th symmetric and the 3rd-10th 
unsymmetric modes which are calculated 
should be printed. 

Logical variable whose value is .TRUE. if 
the boundary conditions should be calculated 
but not stored. 

Logical variable whose value is .TRUE, if 
the planform is symmetric. 

Logical variable whose value Is .TRUE, in 
conversational use. 

Logical variable whose value is .TRUE, for 
conversational use when Input is from the 
terminal. C0NV2 controls conversational 

prompt i ng . 

Unit number for conversational Input. 

Unit number for conversational prompting. 
Unit number for input. 

Unit number for output. 

Unit number for boundary condition file. 
(JJMAX+1)/(MM+1), where MM = the number of 
spanwise control points on the entire wing. 
Number of times that subroutine FLPDWN 

should intergrate to find the flap downwash 
mode or modes. Subroutine FLPDWM has not 
been modified yet to be compatible with the 
boundary condition program, however. 
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OTHER ARGUMENTS: 

(ALFA) Storage space for downwash modes. 

(ALFACS) Storage space for flap downwash mode. 

5.3 SUBROUTINE SLOPES 

This subroutine Is called by subroutine BOUND to 
compute user suppl led modes. 

See sections 2 and 4 for conversational and batch 
use . 

The arrays A, B, C , L, G, E, and H should be 
dimensioned as large as the number of table entries. 
Arrays CHI and ALF should be dimensioned as large as 
the number of table entries plus 1. The dimension of 
the arrays TEST and CHIDUM should be as large as PP 
and 2*PP, respectively. The dimension of ALFA3 must 
be the maximum of 2*PP and NSMAX. NSMAX is the „ 
maximum number of spanwise stations at which data 
should be given. NSMAX is a dimension for ETA, 
ALFA2, and ALFA3 . NPPMAX is the first dimension for 
ALFA2 and the maximum allowable value for PP. 


INPUT VARIABLES: 


U5 

U6 

UCO 

C0NV2 


( ETACP) 

(CHICP) 

MMU 

PP 


Unit number for input. 

Unit number for output. 

Unit number for conversational prompting. 
Logical variable whose value is .TRUE, for 
conversational use and when input is from 
the terminal . 

Spanwise stations at which boundary 
conditions will be determined. 

Chordwise stations at which boundary 
conditions will be determined. 

Number of ETACP values. 

Number of CHICP values. 


OUTPUT VARIABLES: 


(ALFA) Downwash boundary condition at the control 
points defined by ETACP and CHICP. 
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5.4 SUBROUTINE USLOPE 

This subroutine ts one which is to be supplied by the 
user in the case that he wants to define the 
chordwise slope distribution at some particular 
spanwise station by a method unavailable in 
subroutine SLOPES. This subroutine will be called by 
subroutine SLOPES whenever |TYPEC<0 or ITYPEOB, The 
subroutine should return the actual slope values not 
Including the twist. Subroutine SLOPES will apply 
the minus sign (because the induced downwash must 
equal minus the slope) and then add in the twist. 

INPUT VARIABLES: 

Unit number for input. 

Unit number for output. 

Unit number for conversational prompting. 
Logical variable whose value is .TRUE, for 
conversational use and when input is from 
the terminal. C0NV2 controls 

conversational prompting. 

The number of the spanwise section, 

lOKNSPSEC . 

The spanwise location of the station at 
which USLOPE is to calculate the chordwise 
slope distribution. 

Number of (CHICP) values. 

Chordwise stations at which boundary 

conditions will be determined. 

In Subroutine USLOPE, this variable will be 
<0 or >6 upon entry and can be used as a 
parameter If desired in subroutine USLOPE. 

OUTPUT VARIABLE: 

(ALFA3) The value of the slope at the PP stations 
defined by (CHICP) . 


U5 

U6 

UCO 

C0NV2 


N 

ETA . 


PP 

(CHICP) 

ITYPEC 
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6 INPUT FILES 

The following disk files are read by the program. 
The AMES' TSS version of the program issues Its own 
DDEF commands for the files# so none need be given. 
For other systems appropriate control cards will have 
to be supplied for units 8# 9# 11# and 12. 

6.1 GEOMETRY FILE 

This file Is a variable record length file and is 
read from unit 7. 

The first record contains identification and title 
information including number of control points and 
Integration points. 

The next record contains the chordwise control 
points# the array of indices from which the soanwise 
control points are derived# the tangents of the wing 
edge sweep angles at the Integration stations# etc. 
For a complete description of this file see ref. 2. 

On the AMES' TSS system this file has the name 
GEOM.XI where I is the numerical value of ID1. 

6.2 Aerodynamic Influence File (AIM file) 

This file is a variable record length file and is 
read from unit 11. A detailed description is given 
i n ref. 3 . 

The first record contains identification and title 
information plus information about the size of the 
matrix and location of spanwise and chordwise control 
points . 

The second and subsequent records contain the 
influence matrix itself. This file is generated by 
the influence matrix program. 

On the AMES' TSS system this file has the name 
AIM.XI.XJ where I Is the numerical value of ID1 and J 
is the numerical value of 1 02. 

6.3 Identification File 
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This file is read from un i t 9 and rewritten on unit 9 
and contains identification numbers in binary form. 
The third number on this file (IP3) is incremented by 
1 and then the file is rewritten using the 
incremented value of 103. IP3 is the identification 
number for the boundary condition file. 

On the AMES' TSS system this file has the name 
I OF I LE . 
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7 Boundary Condition File (BC file) 

This file is a variable record length file written on 
unit 8. A detailed description of this file is given 
in ref. 2. 

The first record contains identification and title 
information plus information Identifying the type and 
number of symmetric and antisymmetric cases. 

The next NSYM records are right-hand sides (i.e. the 
(BCS) downwash modes calculated by the program) for 
symmetric cases. NSYM equals the number of symmetric 
cases. The next NASYM records are the right-hand 
sides (i.e. the (BCAS) downwash modes calculated by 
the program) for antisymmetric cases. NASYM equals 
the number of antisymmetric cases. In the case of an 
unsymmetrlc wing there will be NSYM + NASYM 
right-hand sides. 

On the AMES' TSS system this file has the name 
BC.XI.XK where I Is the numerical value of 1D1 and K 
is the numerical value of ID3, which is determined 
from 1DFILE at the time the program Is run and is 
found in the program output. 
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8 SAMPLE TERMINAL SESSION 

Given below Is a sample terminal session Illustrating 
the conversational use of the boundary condition 
program on the Ames’ 360/67 TSS computer system. 
This session mfght also be useful for the batch user 
to study. The program used a previously created 
geometry file to obtain the Information needed to 
calculate the boundary conditions and determined 
boundary conditions for 4 symmetric cases and 2 
antisymmetric cases. 

The symmetric cases consisted of (1) uniform 
downwash, (2) pitching about the origin, (3) and (4) 
user-supplied modes. The data supplied for each 
user-supplied mode was identical except that linear 
spanwlse interpolation was requested in the first 
user-supplied case, while C0MM2 spanwlse 
interpolation was requested in the second. The data 
in each of the two user-supplied cases consisted of 
twist and camber distributions at 4 spanwlse 
stations: eta = 0.0, 0.2, 0.6, and 1.0. The twist 

In degrees was determined by the following equation: 

twist = -10*eta 

Values derived from the above equation had to be 
converted to radians. 

The camber distribution at eta » 0.0 was that of an 

NACA, 5-digit, 230 mean line, whose equation is 
z/c = 2. 65 95 *(. 114714984* CHI-. 6075* CH 1**2+ CHI** 3) for 
01CHI<.15 

z/c = 2. 6595*(1.-. 00830377* CHI ) for .lS^CHia. 

This section was defined to the program in terms of 
these polynomials. 

The camber distribution at the remaining spanwlse 
stations were all parabolic arcs given by the 
equation 

z/c = .1*CH!*(1.-CHI ) 

This camber distribution was defined to the program 
In different forms at each of of the remaining 3 
stations. At eta ** 0.2 the polynomial representation 
was used. At eta * 0.6 and eta = 1.0 a table of 11 
values derived from the above equation was used to 
define the section. At eta = 0.6 C0D1M2 

interpolation and differentiation was used while at 
eta = 1,0 a spline fit to the curve (SUBROUTINE 
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SPLIN1) was used and differentiated. 

The antisymmetric cases consisted of the rolling mode 
and a single user-supplied mode. The user-supplied 
mode was simply a case of twist only with linear 
spanwise interpolation between the two stations eta ~ 
0. and eta = 1. The twist in radians was simply 
equal to eta so that this mode was identical to the 
rolling mode with a minus sign. 

The actual terminal session is reproduced below with 
additional comments added In parenthesis. The output 
from this session was directed to a disk file named 
OUTPUT. BC.N1. The ODEF (control card) was created 
automatically by the program and the file name was 
computed using the value of ODISK. The contents of 
this file were printed and are given in appendix 1. 

LOGON user id, password, terminal id 

AMES USYSLIR 

JOBLIBS SYSULIB 

JBLB MEDAN 

CALL RC$ 

(The boundary condition program Is now in control.) 

ENTER BATCH 

F 

ENTER INPUT DEVICE NUMBER 
0 

ENTER ODISK 

1 

OUTPUT (is) ON FILE ... OUTPUT . RC . N1 .. . 

ENTER 101,102 
4 

(Now the program opens the geometry file whose name 
is GE0M.X4 and reads this file.) 

(if an appropriate input device number were entered 
and an appropriate DOEF command given, then the 
remaining prompts would have been suppressed and the 
program would read data from some previously created 
data file.) 

ENTER (BC5)/( RCAS)/PR|NTA,PRINTR,UST r 'RF 

TTFFFTTFFF 

FTFTFFFFFF 

TTF 

(Mow the program prints some heading and other basic 
Information and then enters subroutine BOUND.) 
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ENTER C.M. POS./BREF 

0 . 

(Now BOUND calls subroutine SLOPES*) 

ENTER NSPECS, ITYPES 
4 0 

ENTER ETA, TWIST 

0 . 0 . 

ENTER ITYPEC 
2 

ENTER NPOLS 
2 

ENTER CHIMAX, SCALE, CHIREF/ 

POLYNOMIAL COEFFICIENTS 
0.15 2.6595 0. 

0. .114714984 -.6075 1. 

ENTER CHIMAX, SCALE, CHIREF/ 

POLYNOMIAL COEFFICIENTS 
0 2.6595 

0. -.00830377 

(The above data constitute the equations for the NACA 
230 mean line. Now the second spanwise station will 
be considered.) 

ENTER ETA, TWIST 
0.2 -.0349066 

ENTER ITYPEC 
2 

ENTER NPOLS 
1 

ENTER CHIMAX, SCALE, CHIREF/ 

POLYNOMIAL COEFFICIENTS 

0 .1 

0 . 1.0 - 1.0 

ENTER ETA, TWIST 
0.6 -.1047198 

ENTER ITYPEC 
4 

ENTER TABLE SCALE FACTOR, Al(OR XKC), 

AN( OR EPS), AND THE TABLE. ( CH I =9 9 . STOPS ) 


0 . 

1.0 

0 . 

0 . 

0.1 

.009 

0.2 

.016 

0.3 

.021 

0.4 

.024 

0.5 

.025 

0.6 

.024 

0.7 

.021 
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• :■$' ' 


0,8 .016 

0. 9 .009 

1 . 0 . 

99. 

(The above Is a table of the parabolic arc.) 
ENTER ETA, TWIST 
1. -.174533 

ENTER ITYPEC 
6 

ENTER TABLE SCALE FACTOR, Al( OR XKC), 

AN( OR EPS), AND THE TABLE, ( CH I =99 . STOPS ) 


0 . 

0 . 

0.1 

.009 

0.2 

.016 

0.3 

.021 

0.4 

.024 

0.5 

.025 

0.6 

.024 

0.7 

.021 

0.8 

.016 

0.9 

.009 

1 . 

99. 

0 . 


(Now control leaves SLOPE returning to BOUND and then 
calls SLOPE again for the second user-supplied mode.) 
ENTER NSPSECS, I TYPES 
4 1 

ENTER XK 

1 . • - ' 

ENTER ETA, TWIST 

0 . 0 . 

ENTER ITYPEC 
2 

ENTER NPOLS 
2 

ENTER CHIMAX, SCALE, CHIREF/ 

POLYNOMIAL COEFFICIENTS 
0.15 2.6595 0. 

0. .114714984 -.6075 1. 

ENTER CHIMAX, SCALE, CHIREF/ 

POLYNOMIAL COEFFICIENTS 
0 2.6595 

0. • -.00830377 

ENTER ETA, TWIST ' 

0.2 -.0349066 

ENTER ITYPEC 
2 
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ENTER NPOLS 
1 

ENTER CHI MAX, SCALE, CH I REF/ 

POLYNOMIAL COEFFICIENTS 

0 .1 

0 . 1.0 - 1.0 

ENTER ETA, TWIST 
0.6 -.1047198 

ENTER ITYPEC 
4 

ENTER TABLE SCALE FACTOR, Al( OR XKC) , 

AN( OR EPS), AND THE TABLE. ( CH I =9 9 . STOPS ) 

0 . 1.0 

0 . 0 . 


0.1 

.009 

0.2 

.016 

0.3 

.021 

0.4 

.024 

0.5 

.025 

0.6 

.024 

0.7 

.021 

0.8 

.016 

0.9 

.009 

1 . 

0 . 

99. 


ENTER 

ETA, TWIST 


1. -.174533 

ENTER ITYPEC 
6 

ENTER TABLE SCALE FACTOR, Al(OR XKC), 

AN(OR EPS), AND THE TABLE. ( CHI =99 .STOPS) 

0 . 0 . 

0.1 .009 

0.2 .016 

0.3 .021 

0.4 .024 

0.5 .025 

0.6 .024 

0.7 .021 

0.8 .016 

0. 9 .009 

1 . 0 . 

99. 

(Now control leaves SLOPE returning to BOUND and then 
calls SLOPE again for the user-supplied antisymmetric, 
mode . ) 

ENTER NSPSECS, ITYPES 
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2 

ENTER ETA, TWIST 

0 . 0 . 

, ENTER ITYPEC ■-> 

0 .. .. 

ENTER ETA, TWIST 

1 . 1 . 

ENTER ITYPEC 
0 

(Now control returns to BOUND and then to the main 
program.) 

ENTER INPUT DEVICE NUMBER 
-1 

TERMINATED: STOP 777 

(The operating system Is now in control.) 

PRINT OUTPUT. BC.N2,PRTSP=EDIT,STATION»RMT05 
PRINT BSN*????, ??? LINES 

LOGOFF 
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OUTPUT FROM SAMPLE TERMINAL SESSION 



DETERMINATION OF. Thin LIFTING S UR F A-CEL -B .QiLiiD A R. T.~CCiiDXT J-0 NS.. 


AjJ 

o 


TITLE a RECTANGULAR M«G 
ID1 = A 

102 ... s . 0 

103 a 7 

RR 4 T = 1,00000 

CBAPRR a 1,00000 

rHOP0*ISE CONTROL POINT locations 
0,000000 
0 , \UbUii7 
0 .R 0 O 000 
0, 6^3553 
1,000000 

SP a nwJ 5 E CONTROL POINT LOCATIONS 

<*■ , 

0 , Pbb 0 23 
0,707107 
0 .S 000 OO 
0 ,25*1? 1« 

- 0,000000 


AR = 2 13-13-75 


(>Cs 

T T 


E F T f FEE 


HC AS 

(T T F T E F F F E F 

P R I N T 6 , PKTfvT*, ft STO^E 
T T v 



UNIFORM OCJWMWASH 


i . 0 . 0.00 ooo 

r. ooooooo 

r. 000000,0 

l'. ooooooo 

1,00000 0 0 

i. ooooooo 

1 '-0000 0 00 

!'. OOOOOOO 

i\ ooooooo 

1 ,0000000 

1 ,0000000 

- 1 .00 0 0000 - 

J-t-i 00-0 0 0 0 0 

i ’.ooooooo . 

1 ,0000000 

1 ,0000000 

1 *. ooooooo 

t 1 '.ooooooo 

i . ooooooo 

• 1 , OOOOOOO 

1 ,0000000 

• r'. ooooooo 

1 '. ooooooo 

l \ ooooooo 

-• .. 1 ,-On) n o ft 0 0 

V. 0 0 00 0 0 0 

t'. ooooooo 

. 4 '..00 0000 0 

: 1 ‘.OOOOOOO . 

1 , , 0000.000 

PITTING MODE 

i. 


r ' * 


C , M , POS’./.BREF 

s 0,000000 




0,0000000 

o'. J a6aa66 . 

0-.5000000 

. O'. 8 5.555.i3 -- 

1-, ooooooo 

0,0000000 

o'. 

O', <5000000 

0 ! , 85355^3 

1 .OOOOOOO 

.. ^ « 

0,0000000 

o'.U b^U 

0, C$000000 

0.8535533 

1 .ooooooo 

0.0000000 

o'.l«o«a fc6 

O'.sOOOOOO 

- 0*. 8535533 

.1,-0000000 

0,0000000 

O'. 14^4466 

O'. 5000000 

0.8535533 

1,0000000 

0.0000000 

••• 

o'. 5000000 

0.8535533 

1 .ooooooo 



u sfr -supplie d_s y_m h ejr. ic ... mode 


2329.6 1« _ . 

■ 0,2400104. . 

.0,1685513 

-0,0900664 

-0.0743155 

*fl ,2238 i 79 

-0,2223796 

-0,1511248 

-0,0005781 

-0.0567508 

-0,2092727 

.0,1943331 

.0.1234033 

mQ, 0527584 

-0,0283894 

_ _.-0, 1.83320? .. .... 

■3,1579772 

. -0.0872664 

-0,0165554 

0.0087901 

• 0. 144400a 

•0,1158831 

-0,0451724 

0,0255383 

0,0540544 

-0,3050843 

-0.00298U4 

0,0220839 

0,0220839 

0,0220839 

USFR. SUPPLIED S 7 MMF TR Ir MODE 




. ji 0,234 0535 

-0,2399799 

• 0,16.85526 .... 

•0,0900581 

.0,0745003 

-0,2269311 

-0,2222926 

• 0, 151 12a9 

-0,0805547 

-0,0572790 

-0,2120129 

-0,1942565 

-0,1234068 

-0,0527377 

-0.0292703 

.■0,1.85 Oil 6!. 

-0,1579284 

-0,0872711 

-0,0165292 

0,0084959 

-0.1379284 

-0,1161122 

-0,0451886 

0*, 025627 0 

0.0553199 

-0,3050843 

• 0, 0 029644 

0‘, 0220839 

0 *. 0220839 

0,0220839 



RQLLIN£_.MQa£ 


0.^659258 

O', 9659258 

O', 965-925.8 

Q ’.9659258 

0.9659258 

0,8660254 

0*. 866 0254 

0,8660254 

0 t .66fe0254 

0,8660254 

0.7071068 

O'. 707 1068 

0,7071068 

0*. 707 1 068 

0,7071068 

n 50 0 0 0 0 <L_ 

0* 5000000 

0^50.00000 ........ 

.... . 0.',.500 0 0 00 

0.5000000 

0.2508190 

0*. H 5 00 190 

O‘.250019O 

O', 2588190 

0,2588190 






USER-SUPPLIED A NT I 

•SYMMETRIC mode 



-0, 9659258 


-0,9659258 

-0,9659258 

-0,9659258 

•0,8660254 

-0.8660254 

-0,8660254 

.0,8660254 

.0,8660254 

•0,7071068 

-0,7071068 

*0,7071068 

-0,7071068 

-0,7071068 

•0,5000000 

-0,5000000 

'-0,5000000 

-0.5000000 

-0,5000000 

-0,2506190 

-0,2580190 

• 0 , 2508190 

-0,2508190 

-0,2588190 


MODES STOHED AND FILE CLOSED 
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c • • 
c.. 

' - — 

c.. 

c.. 


...SURFACE PROGRAM, 
...AND. OPT I ON ALL V . 


\ j 4 

\n 


— A-XHlAL. W iNG-tl F-T-I NS. 

I T USES INFORMATION STORED ON THE GEOMfTRY FlLF 
. ON THE INFLUENCE MaTRIX FILE, ThE BOUNDARY 

. »C 0 N D.J J JO MS. .. J-H A T — IX— C-A L C.U L A T£S ARE ST Qfi£0_0lN A ...8.0 UND. A RY.._C 0 N D I T I Q ji F ILE 

...AN EQUATION SOLVING PROGRAM WILL REaD THE INFLUENCE MaTRIX AND 
.BOUNDARY CONDITION FILES AND DETERMINE THE COEFFICIENTS IN 

1 h£l,EX~P A N S 1 0 N FOR nElTA-r.P, , _ 

.DOUBLE PRECISION £Ta, SThETa 
INTEGER UCI.UCO#U?.U5R,U6,U8,UJl,U2O,O0I«K f U7 

I N.T.E.G' £R-_I.LI LEL, P, » P P .^UTiLS 1T_M ...... _ , 

- LOGICAL C0NV,B.aTCH,L0G1,CQnv2,BcS,BCaS # p"rin.Ta,pRINTB 
- lOGIC'AI STORE, NS TORE * F| AP, SYMM 

REjAJ :tyXJi^LAMQA£__ • ■ ■ • . ■ ■ _ 

DIMENSION TITLE- (26), CHICP ( 20 ) , NI NpEX ( U.7 ) , ET ACP ( M7 ) : 

.. DIMENSION. FT A (383) ,STHFTAf?83)-*TANLFL(i83) »TANLeR(383) 

- DIMENSIO N TAN TE L f 3 ,8 3.) ,.TANT ER ( 38.3UCQRDIP (393 ) « XSILI F.( 383) ' 

DIMENSION cHIFPI(381)»C2IP(383) •BCS(iO)#8Ca8(iO) 

DIMENSION, X0.UM2C2) .XDUMb(b) ,NDUm9(R) , NDUM2(2) * 

EQUIVALENCE (X0UM2, XdUM6,NDUM9.NDUM2) 


Data ■.UCj*UCO,U5,U5R,U6,U8,uil,U20/5,6, 
data* coNV 2 / t false,/ 


5, 5, 6 , 8, li,9./,U7/7/ 


C* . • • ' - . • .... • • ! ' 

C UCI s UNIT NUMBER' FOR CONVERSATIONAL INPUT* 

IL^^Ji£P^.s UNIT. NUM BER FOR CON VERSATIO NAL OUT P g T , _ 

c us s unit number for input of tmE first card and’ 

C..... • THE. .UNIT number FOR INPUT UNDER CONVERSATIONAL USE. 

. C.« ... . . U5R S,„ UNtT n umber Fqr tNPUT (EX CEPT F OR TmE FjRST _£aRD) 

C , UNDER BATCH CONTROL, 

* T ME UNIT NUMBER FOR REGULaR OUTPUT UNDER EITHER 

C. .... -9-A.TO OR_CQNVERSaTIONaL contro l . 

C up.. B UNIT NUMBER FOR BOUNDARY CONDITION FILE. 

C U 7 , = UNIT ; number for reading the geometry file 

*_aili s . UMT NU MBER FOR READING TH E INFl U SJi C £_ » Al ft I X F ILE 

C U20 , 5 UNIT NUMBER FOR IDENTIFICATION NUMBER FI|E. 

C - • 

— t .... . . the FOLLOW I NJ1_IS_NF,CL§^A R T EM -IS.S 



c 

US 5 5 

c - 

WRITE (UCOf-feO 04 ) 
c..,..oh batch use'. 

c . .. . 

-C THE FOLLOWING ALLOWS The. USER To CHOOSE CONVERSATIONAL 

HEAOfUSf 5000) BATCH 
- C0.NY-s.NOU9A.TCH 

c 

20 CONTINUE 

— — IF.f-BA.TCfi).. GO TO UQ . 

w R T T E f U C 0 » 6 C 0 1 ) 

READfUC USOOi) 05 

IF. . f.US._ * E Q ... 04 .US a. DC I . 

C On V 2=U5 , EO . iJt I 
C 

-X. ... ... U5.LT..Q IS ..FOR . TERMINATING EXECUTION CONVEHSA T IONaLL* 

C 

lFfU5.LT.0j STOP 777 

- C^,.,.IF (CPNV 2 ), THEN ALL The. input w I L L a£ FROM The terminal 

30 CONTINUE ' 

WRITE (UCO/6005) 

c 

C ODISK IS FOR The a*ESi TSS VERSION of THE PROGRaH 

C..., .ODISK MUST fit AN INTEGER FROM n T 0 9. 

— C.» ti j.* iF .QPl.S.K s. Of OUPUT WjLL BE On.UNjT b t jF OPiSK jS BETWEEN 0 AND 1 A 
C, # .,.THEn OUPUT W t LI St ON UNtT a AND ON A F j L E NAMED OUTPUT, BC. NX 

C WHERE X s VALUE OF ODISK', THE 4 M£S i P»OGHamS OBEY AND CVHT aRf 

C.. .USED, T0..0I.N|RaTE AND. GIVE THE NEEDED DOEF and RELEASE CUMManqS TO 

C TSS OPERATING SVSTEm’. 

C 

R F.A P. .( UC I f 5 0 0 .1 ) . ,0 D,I S K 

U6 = 6 

IF COpISK f E3, 0) GO TO 50 

Up s._ ji 

ODISK - MOOfMAvOfl .ODISK) > 1 0 ) 

WRITE fUCO,600H) OpI SK 

CALL 0 S tV(lfc,J6HRfcLEAS£ FTOaFOOl ) 



CALL CVftT(ODlSK,l, 

1 tar(i|)OEF FTOaFOOl > » OUTPUT. BC.h' , tl #6X) 

- ? ALFA $ R, 5 

Cf-ll 03fe V( I'd . A IF O • • • . ' - 

REMIND U* 

. 50. TO Sfl ' 

«f! ro^l I «Ut > ....... 

USsUS* 

50* C 0 ‘V T I s U £ : ..... .... • ' 

T f t r 0 N V ) w R I T E f, U C 0 1 b 0 0 0 ) ■ ' 

IF f r.v -vV) Rfc AO ruC i 1 500 1 ) -X0Dl f IPO? 

IF cBaTC^KEaP-CUS.SQOI } '..ID-01 . IDD2 .. . 

• IFf 1001 .IT-.OISTOP 

c • 

C . * , ... I s - ! ..THE PqP- 1 0 .. VERSION. IDD I AND IDD2 ARE • USED TO DETERMINE THE 

r p.FOmFT »V ANiD INFLUENCE MATRIX FILE NAMES', SAME -FOR AMES! TSS VER'. 

C «'«»«« T N AM Y Case, Twt VALUES R.E4D in WILL, if GREATER Than zero, 
c ..... «e checked aga i M s.r'.T.H.t, values on. the- ..geometry and . influence. 

C • • . . . M A T R I X FILES AMD A STOP OR PAUSE WILL Ft EXECUTED IF THERE 
C IS NO agreement, 

CV.....-IE IDD.2 ,.LK. 0, THE... INFLUENCE _m.aTRIX • FILE WILL NOT rE USED’. 

c . • . : 

r .' ■ 

C..,.,,.F0R AMESi JSS VERSION,: ISSUES. pQEF COMMANDS TO OPERATING SYSTEM, 

r : . 

' Call GEMFIL(TDDI) • .. .. • . 

REaDlU 7 . ) 10 J ,RP r MM 4 NDUM2,UNSYM f NOUMg.MREF, JJH a* , NFL A PS, 

1 T I T l_ f f N T I T L 

r . •'.... ‘ " ' 

C,.,,,,PP • _jt_NUMRER ■ OF , CHORQNISE CONTROL POINTS 

C, .Vm^' ’ s NUMBER .'OF SPANRT.SE CONTROL POINTS', 

C UnSYk s NON-ZERO for UnSvmmETSJCAL wiN G ', • • 

o, jlj J*MREF „_s..REf.ERENC.E_-.TNTEGE*R FOR DETERMINING THE SPaNmJSE 
P,,.,. CONTROL 'POINTS FROM ( N j N.OE Xj AND (ETA), ThESE 

CONTROL, POINTS WILL rE STORED IN (ETaCP), " 

C.....TJHAX = JT_b EJ. N U M 0 E R 'Op SUNWISE INTEGRATION POINTS', THE 

c spanwtse control points are a subset of The • ‘ 

c , . . sPA'Nwise .integration points'. l- 

C * • *.*.«..N.F.L,AP-$._s N.U.Mg.E.R„...0 F _.C.O N T RO.L- SURFACES fOtl* OR 2). 



C,.,.. (TITLE):: ARRAY OF ALPHANUMERIC TITLING INFORMATION, 
c NTlTlt s ThE NllMRfc R OF WORDS IN ThE TITLE , 

C LOG I s Id ,NE. IDOl .AN0, IDD1 ,GT, 0 

IFfCONV, aNO.LOGI JPaUSE * 10 NUMBER CONFLICT FOR GEOM FILF » 

IF cBATCh.ANO.. LOG 11.. C.ALL STQP2(U6,» JO NUMBER CONFLICT FOR 
1GF.0* FILF i t FLOAT ( IDl ) ) 

w £ A D C U 7 ) fCHlCPfl) t T = 1. #PP) , {NInOFX(I) . I«lr" M > . ( T ANLEL ( I ) , Is \ , 
_ UJMA-X.l*JTANLER.airl=li JJMAX)^ (TANTFLCI) * Isl, JJMAX) , 

?(TANTER(I J , Isl , JJMAX) , f ETA cl) , 1*1 , JJMAX) , 

3 f ST MET AC I) , Isj # JJMAX) r CXSILIRCI)*Isi»J'I h AXi, 

u CCORD.IPI n.«.Is±*J.JJ)AX.) ,B-RAT.CBARBR 


fCHlCP) s ARRAY OF CHORDWlSE CONTROL POINTS SUCH THAT 

C 0.0 .LE. CHICPCI) ,LE. I. 000 . 

KU'OfcX) = Integer ARRAY from which THE SPARSE CONTROL 

POINTS ft RE Ofc«mO FROM THE SPaNWISE 

c INT£GR aTION .points. 

MAN) El). CTaNiER). cTanTEL). AnO f TANTER) ARfc ThE 

c lt tangents of the wing edge sweep angles at the 

C..^., .JLNlEGRtUQN stations. 

(ETA) = T h E SPAmWjSE INTErRaTIOn STATIONS NON-0 i HENS i ON AL- 

39 r I ZE D R Y the EFFECTIVE SEMI-SPaN, B2 ((ETa) IS 

C..... . OQURLE PRECISION). 

r csthETa) = c srrt n . *f t * * *? ) ) (double precision), 

c CXSILIP) s LEADING EDGE LOCATIONS AT THE INTEGRATION STATIONS 

__ normal.jzeq_.sv b2, 

C CCOWOIP) s STREAMWISE CHHRO LENGTHS AT THE INTEGRATION STa 

r...., t ions normalized by ep. 

C.....HRAT _ = LATERAL REFERENCE LENGTH (USUALLY ThE TRUE SEMI- 

f Ml( , SPaN)/B2. USUALLV brat WOULD RE 1,0 IF THE WING HAS 

r,,,., NOT BEEN Yawed. 

Cij •» .CBAMPR s_ LONGITUDINAL REFERENCE LENGTH (USUALLY THE MEAN 

> . , , . . GEOMETR IC ChORD) /LATERAL REFERENCE LENGTH. 

r. 

IFfNFLAPS.NE.O) READ ( U7 ) XDUM,LAM0AC.X0UMa,ETAl,E.TA2 # . 


jXrUMb. CrHlFFl(I ) t 1=1 . JJMAX) , CC2l p (I)» I« 1 »JJHAX> 

c 

C.^^iLNE INFLUENCE MATRIX FILE WILL NOT RE READ IF IDD2.LT.0, 



c 

IFf IDD2.LE.0) GO TO HQ 
C 

r • . 

... C . ... . .*.F OP. a*£S» TSS VERSION,. . ISSUES ,DD£f COMMANDS. 

,c ‘ ‘ ■ 

C&LL, AI^PILCIDOI, 1002) 

RE. A D ..( U 1 n . I D1 1 1 0.2 ..MDUM , T I.TLE-, Nil XL. P.P.,.NDUH # Mri . MRXF 

J NDOMQ, MACH, tPSr XOUM, ( C H I CP ( 1) , I c 1 , p P ) * f N I MOE X ( I ') > I = 1 t MM ) 

INFORMATION FROM the INFLUENCE MATRIX FILE OVERLAYS INFOR- 

c,.., ..Nation from the geometry., file,.; - : . 

LOG 1 s I001 t Ne.O...',AND, IDDi.Nt.IOl 
LOGlaLOGI .OR. f I D02 , NE , 0 .and, IDDa.NE.lD2j 

„• i f t.c o n v ..and *l c,g.u .. p. ajj $ t i o.: num&.eb_xj3m uci ..e.or, influence 

1 MATRIX file ' 

IP fBATCH, ANO.tOGj ) Call 5T0P2(U«,,t ID NUMBER CONFLICT FOR 

llNFLUENC. E M AT. R.IX JLI L E .. L.. • # F l 0 a T. ( 1 0 ?.) ) 

fcO CONTINUE : 

•c . 

i.FrcnNV2i..wKnexUcOfj»jo.o.7) _ . 

^ RE AO f U5 * SO'OClT PCs • 

° Rj AD f U5# 50 00 J 3C AS 

READ (U5» 5000 J PRINT A j RRINTR, mstqRE 

PRINTB-PRInTb.OR.PRInU . , ' 

C " ' 

C.....EACH ELEMENT OFrRCS A'^O 0 CaS CORRESPONDS TO a DIFFERENT DQwnwaSH 
C .CONDITION, A TRUE V a L LIE MEANS TH A T THE CORRESPONDING bOUnDaRY 

c ., ,. .condition rill ue calculated or read from cards and, if 

c nstore is .false... then .written on a file for use in the 

- c ECUaTION SOLVING PROGRAM. both arrays are read and acted upon if 

C » « • • • ThE WINR.-IS SYMMETRIC OR UnSYmmETRic, in ThE LATTER case T m E 
Ci^jjsBOUNDaRY CONDITION will. BE CALCULATED FOR the ENTIRE. WING. 
C....;fhF CASES CORRESPONDING to Each OF THE ELEMENTS aRE aS FOLLOWS — 
r. . 

, c gcsiij . untfo.Rm mode *... 

c RCS(?) PITCHING mode 

r ALMfX8I,ET»)afXSI/f j RAT-XSlCH.)/CHARBR' . 

c WHERE XSICMsC.ENTFR of MASS OR. other REFERENCE POSITION.' 



nonnnbnnon n’n on'hnci^onj^nnpnn 


c 

r 

„.c 

c 

c 

-C. 

c 

C R C ? ( 3 > 

C BCStu) 

C 

.._C_b.CS c 5) — 


OF 1. 

TO THf 
COORDINATE. 


■XSIC M IS TO p£ GIVEN IN The COORDINATE system fixed 
TO Thh; WIng WHILE X S I IS I* Th€ WIND-CENTERED 

COORDINATE SYSTEM,. , ... 

THIS NODE CORRESPONDS TO a N0N-DIMEN9 ION aL PITCH HftTE 
thp solution for This mode gives the ring contribution 
Q STABILITY DERIVATIVES (IN THE WIND CENTERED 
SYSTEM} , 

linear, symmetric twist mode 

_ ALE A ( XSI « ETA) sABSfET A) /BR aT 

parabolic, SYMMETRIC TWIST MODE 

ALFA(XSIiETA3B(ETA/BflAT)**2 

RESIDUAL ft AP DOWN.WaS.H. mode. IF.. THERE aRE 2 flaps, 

THIS MQnE is FOR The symmetric DEFLECTION of both, 

REFER TO *AS A TN DW251 AND TC SUBROUTINE FLPDWN FOR 
FURTHER DOCUMENTATION.. ..the flap DOwNwaSM modes are. the. only 


.T.HE.N. 


BCS(fe)- 

BCS(IO) 


ONLY Ones AFFECTED 9 Y THE MaCH NUMBER, 

TwE value OF 8CS(B) ANO BCAS(S) 'WILL BE IGNORED IF TwERE 

15 NO.. FLAP DATA AVAILABLE .FROM.. THE..GEOMETR.Y FILE,-- - - 

THESE ARE MODES Fqr which Tables ANO/qR certain TYPES 
OF COEFFICIENTS will BE read BY the PROGRAM, 
the MODES are .DEFINED BY ..SIREaMHISE. DISTRIBUTIONS ... 
(COEFFICIENTS OR TaBLES) aT SPaNWISE STaTIONR SELECTED 

pv the user, these modes are assumed symmetric only 

IF THE.wlNG. is. S.Y iM.ML.TRI C SEE_ SUB ROUTINE BOUND FOR 

FURTHER DOCUMENTATION, 


BCAS(l) a NT I •SYMMETRIC UNIFORM MODE (ZERO aT CENTER) , 

BC A S ( 2 ) ROLLING MODE (RIGHT TIP DO^Nl 

_AlFA(XSI.,EtA).sEIA/BRAT. - 

THIS mode CORRESPONDS to a NON-DIMENSlONAt ROLL Rate 
OF 1 . 0 . THE SOLUTION FOR This MQ D E GIVES THE WING CON- 
TRIBUTION . .TO ThE. QUASI-STaTiC, P STABILITY deriy- .. __ 

■ aTIVES (IN THE WIND CENTERED COORDINATE SySTEMi, 

BC a? ( i 5 ANTI-SYMMETRIC RESIDUAL FLAP MODE, IF THE WING IS 

UNSYMhftRIr, There is mo. NEE(5-TCLINV0KE .This. «OoE., 

BC f S ( u ) “ FOR a SYMMETRIC PLaNFORm THESE MODES a«E THE ANTISYMMETRIC 
bc a s c i o ) counterpart of those referred to by bcs(6)-bcsc 103 , 

FgR an UNSYMMETRIC.. WING the treatment OF .THESE MODES 15. 



r> n n n n ^ n 


c THE SAME AS THOSE CORRESPONDING TO BC S ( 6 ) *6CS ( 1 0 ) , 

C 

r PRINT* T RuF- HFaNS that aUU THE DOWNWaSH MODES WHICH 

; ; ; • ___ : ie£. C aLCUL A lED_-JULi-.RE-.-PR mtED. 4 R£GaRDLESS-.OE_J.HE -V.AlUE. Of 

PRINT ai , 

;;:;!p«intb '.true, means that the 5TH-iOTM symmetric and the 

f 3RD- 1 0-T H... U N S..Y. MM E..T.R I-C ...MODE S- JiH.J-C H_.AR.E-- C AL-Cil L A-T.ED~-W.jUC - 

!!!!! be printed'. 

NSTQRE TRUF, means that the boundary conditions will re 
; ; \ \ \ C A LC.UUT E D BUI. AtO t._SI.QRE D . :■ - 

RCSCOsBCSfS) .AND. NFLAPS.GT.O 

B C A. s C i i * B.L.AS f 3 ) . AND ,.nF l. A P S . G.t,.Q. - — -• 

FlAPsBCS(S) ,0R,8CAS(3) 

IF(C0NV2.AND,FLAP) WRITE (UCO, 6021) 

IF (F.LAP) READ f U5» 5Q0l.)_ Li - 

C 

STORFa. NOT.NSTORE 

1.0 3 5 A . . . — - 

■p- SYMMsUNSVM.fcQ.O 

1-1 MMPsMM 

.MM.PAs.MM - - - 

IF f SYMMj MMPs f MM+ 1 ) /2 
IF (SYMM)MMPa=MM/2 

C;;...(ETACP)sLOCATIONS of spanwise control STATIONS, 

JRATIOsC JJMAX+ j )/(MH eF+i) 

DO 85 M*l#M Mp - - 

INDEX s N INDEX fH ) *jr at ro 
ETACP(M)*ETAf INDEX) 

85 ' CONTINUE... 

IFCNSTOrE) go TO ro 

r DETERMINING I.pENT I F I C.A T I ONi NUMBER OF THIS RUN 

'"’’call obey (?2»2ljhddef ftqrfooi # »idfile ) 
rewind uso 

RE.aO(U2C) I.DAi IOB, IPJi TOD 



n n mm td n n n 


103 * 103+1 

RF^lNiD U20 

WRITF.(U2o) IDA, IOB, ID3# IDO 
END file U20 

c FOR A^ESi TSS VERSION ONLY 

.CALL C8EYn6#l6MRtLEASE FT09F001 ) 

,.,,»N5YH - NUMBER OF SYMMETRICAL CASES^ 

.....NASYM = NUMBER OF ANTI-SYMMETRICAL cases. 

N s Y w s o 
N A S v M s 0 
00 B7 1 = 1, 10 
IF (8CSf I ) JF'SYMbNsym + I 
IF(3CASCI) )MASY«sNASym+1 
7 CONTINUE 

.... .writing ThE introductory record of the boundary condition file'. 


FOR AMFSt TSS VERSION ONLY, pCFIL ISSUES 00F.F COMMANDS TO THE 

.....OPERATING system. 

call bcfil noi , ID3) 

►RITE t US) HJirlD3,TlTLE,UNSYM,NSYM # NASYM # BCS,BCAS,PP,CWTYPE, 
1 MM,MMP # MMFa,MREF»SWTVPE, { CH I CP < P) , Pal ,PP ) , f NI NDE* f I) , Isl , MMP) , 
2(ETACPf*) . Ms 1 ,MMP) 

90 .... CONTINUE 
C 

r 

WRITeCUH, 6003) 

WRITE (U6,6r 08) (TITLE(N) .N-jjNTjTL) 
iFriDS.FG.O) *RIT(.'fUfc»6009) 101,102 
IF ( IP3 « NF , 0 ). WRITE <06,80205 I01,ID2,ID3 
NRITf. (U6,feOl0) BRaT,CBaRBR 

IPCFI.AP) iNRITp. fU6, 6022) LL 
wRlTFCUfc,6ftll) 

WRITE (U6,fcr 12) t C H I C p C P 1 , p a i , P P 1 
WRITE(U6,6013) 

WR I TF.(U^ i SO 1 2 J CET 4 Gp(I),I = i, N "M p3 



WRITE (life, 6014) 

WRITE (U6, 6*1 5) PCS 

WRITF.CU6.6016) - - 

WHITE (U6, 61)15) BC As 
WRITS (U6, 60 17) 

WRITFCUfc*feftlS5 RRINTa,PR.INT.B,.NSTORE 
WRITS (U6, 6023) 

Cali pound r pp,chicp#mmp,mpp4,etacp,cbarbr,braT, 

1BCS,BCaS,FLaP, PRINT a , . PRINTS, STORE,. STM*, rQNV, C 0NV2,. 

pud. uco. us» U6» ue, 

3XSiLiP,CDROJP, MIWPEX, TaNUeL, TANLfR, TANTfU# TARTAR, 

AiET A , SThElA, LAMDAC, ETAl, ETA2, CHIFPI, C2I p , 

5JJW-AX, JR A T 1 0 , A t„ F A r LL * A L F A C S ) 

IF (UB.fQ.UCOI GO' TO 120 
end FILE Ua 
120 CONTINUE. 

GO Tn 20 

r. 

c * , * » , input FORMATS 
r 

5000 FORMAT f 6 Q u 1 5 
•p- 5001 FORMAT f j 6 1. 5 1 
>J4 5002 FORMAT f6F 10,0) 

r 

OUTPUT FORMATS 

r 

6000 FORMAT (t ENTER 101,102' ) 

600 1 FORMAT (/////• ENTER INPUT DEVICE NUMBER' ) 

6003 fORMATOHI/ 

1 r cOETERRINaTIOn OF THIN LIFTING. SURFACE.. BOUNDARY. . CONDITIONS » / 

6004 ? formatc' Enter b^tch* ) 

.6005 FORMAT (1 ENTER CJQ ZSK • ) . - 

6006 FORMAT!' OUTPUT ON FJLE , , , OUTPUT ,flC.N * r 1 1 f 3H , « , ) 

6007 FORMAT! 1 fNTfR ( 8C S ) / ( BC A S ) /PR I wT A . PR I NTB « NSTORt i ) 

r , THF FOLLOWING FCRVaT DEPENDS gw t he integer word length in 
r! ! ! ! "characters'. CHANGE aS^REQUIREO for THE COMPUTER BEING USED, 

6006 FQWM A T ( // 1 2h TITLE = ,t 20A43 - 



6009 


60 1 0 

6011 

6012 

6013 

6019 

6015 

6016 
6017 

6020 

6021 

6022 

6023 


C 

C **** 

c 


c 

c 


1 00 
H* 

120 

20ft 

210 

22* 


F0RMATf4H 10!., 6 Y , 1H = , IS/ 

1 9H 102, 6X, ins , 15) 

FORMAT (5H a T , 5X, iRis,Fll«5/ 

1 7b. CBARbo, i'i, 1 He , F11.5) 

FORM A T f/ *5H ChORO^ISE CONTROL POINT LOCATIONS ) 
FORMAT (Fp3*6) 

FORMAT f/35u SPANMSE CONTROL POINT LCCaTtONS ) 

FORMAT ( / m h BCS 3 

FCRm*T(1X»10(L1MX)) 

FORMAT ( /5H BCAS ) ' 

F0RmaT(/P9H PRINTa, prjntb, n s t o r E ) 

FORMAT f4H IDtf fcX# 1H = , I5/9H If)2, 6X, IHe, 15/ 

19m ID3, by, 1 ms, 15 ) 

FORMA? ( i ENTER LLI/1 
F0RMA.T(3M LL, 7X, t H* t 151 

FORMAT (■/////) 

E NO 

SUPRQUTINE CODI«ifyi,Xl,wi # T, ANS.NA, XK) 

A CONTROLLED DEVIATION I TERPOL AT I ON METHOD 


DIMENSION X I ( N I ) , Y T f N I ) , T (N A 3 » A N S f N a ) 


N-N I 

SIGN = 1,0 

If (XI(NJ) .LT.Xlf 1)3 SIGN s .1,0 
OP 910 IE=1,NA 
XsT( IE) 

I«i AN-g) i i 0, 120,200 
vs y y f n 3 
GO TO 900 

V s CYI(2)-YI<1))/rxi(2).XTCi))* CX-XHl)) + V I r t > 
r, 0 T 0 9 0 0 
•I 3 1 

Ir ( SIGN* ( Y I r J ) - *n 230 f 220# 250 
v' e v j ( J 3 
GH TO 900 



S*I 


230 J a J + l 

lFfJ-Nl2l0,2)0,'250 
25Q lFfJ*-2) 120, 155,260 
155 ,T s3 
JJ = 1 
GO TO ?8S 

' 260 lF(J-M)2M0,26b,270 

265 J a N*1 

• JJ.. a 2 

GO TO ?65 

270 Vs (YI CN)-YI rN-1 ) ) /CXI (N) -X! <N-l) ) * f X-X T ( W- 1) ) ♦ Y j r N- 1 ) 

GO TO 900. . 

280 JJ a 3 

?65 If fN«-3)290, 290,2«5 

290- J = 3 - ; 

295 5 s J»1 
M a K-l 

L s J*t ' 

Ai s X-XT(M) 

A 2 = 'X-XJt(K) 

--.S3- 5 X-X! CJ) 

AL = (X-XI (•<) )/{XI f J)-XI(KJ) 

S a AL*Yl(J)+(j,0-ALJ«Yl(K) 

C 1 s _ A3* A2/ ( ( X T ( 01 *X J C K) ) * CX I ( H) -X J f J ) ) I 

C2s Al*A3/C(XlfK)-Xl(M))*-(XI<K)-xi(J))) 

C3# A2*Aj/( CXICJ)*XI(P))*(XK J)-XJ CK) )) 

P I * C I * Y l.f.H5 +£2 *Yi CKi +C3* Y J ( J ) 

IF fN-3)305,305> 310. 

3f)5 P2 = Pt. ' 

GQ-.m - 3 t 5 

310 A« a X-X! Cl. ) 

C« = A«*A3/ r C(XlfK)-XHJ))*(XI(K)-XTfL))5 

CSe A 2* A4/H XI ( J] -X 1 £« 3) * (XI (J)- XI CD) 3 

C 6 a A3*A?/C(XICU-XlfKn*CXIfL)-XTfJ))) 

P2 s r«*YT fK)+C5*YT ( J)+C6*YI (L) 

3i5 GO TO .0320, 330, 350), JJ 
320 °2 = Pi 

AL s CX-XTt ln/(xlC2)-XICn) 

5. « AL* V I C 2 ) ♦ 1 1 • 0 * A L ) * Y I f 1 ) 



Pis S + Xk*(P2-S) 

GO TO 550 

33c PI s p 2 

Al s fx-xj (r,-i ) )/(VK>o.V! (N-n ) 

$ S A U * Y I ( M +( 1 ,C-aL j*vi (N-l) 

F2 s S+ xK*C p l-S) 

35o FI s A8S(P1.S) 

F2 s AfS(02-S) 

IF fEl+E2)400,4G0,4l0 

Mfl r V s $ 

GO TO 000 

41“ p T s fM*At)/(El*At+(l,0-AU*e?l 
Y s 6T*P2+ (1 ,0-BT)*Pl 
O0 o A^SflElsV 
91-2 CONTINUE 
9FTURW 

£vr> 

SL'eHOUTiNfc. i»QU*D c PR #CHlcP,MMP # ^MPA # fcTAC° # C6ARBR f BRAT r 

j 9 c s , p c a s » f i a p » print a » points, store, symh, honv, conv?, 

2 < 1 C T , uro, U5 . t'fe, 1'8 , 

JXSIUIP.cOKpIP, f»I^DEX. TaMLEL, TALLER, TaNTeL, Ta*T£R, 
mFTa, STheTa, LAKDaC, ETAi, ET a 2 , CHIFPI, C? I p , 

SJJHA x, JRATJO, Alfa, h , alfACSi 
OOU pLE PRECISION E T a , STHETA 
IwTFCER P, PP, UCI, UCO# 05, Ufe , U8 

L C c I C A L BC8, 8C A S , PRinTa. R»IM8, STORf, CCNV. CONV? 

LOGICAL FUR, SYMM, FLPCaL 

DIMENSION C H IC P (RR ) t ETACPfMMPJ# N I N[)E * ( M pP ) , pCS(tO), BCASflO) 
DIMENSION A|FAfPP, MMP), A| FaCS(LL» pp ' *R p 5 

DIMENSION X.SILIP(JJ.MaX), CORDIP ( JJM A X j , T A^LEL f JJ M AX ) , 

1'TanlER? JJNaX), TaNTEK JJ«aX) , T a^TER ( J JM A X) , ETa(JJKaX)# 

2STHFTA f JJNAX) , ChIfPK JJMaX) , C2IP(J.JMAX) 



FlPCAL.8 tkue.. 




10 

r OnT I MJF 






M M !J 







IF 

( 5CS( 

1)3 

GO 

Tn 

.10 0 

2? 

IF 

( «CSf 

?) ) 

GO 

TO 

200 

23 

IF 

( PCS( 

3)) 

GO 

TO 

300 

24 

IF 

( «CS( 

«)) 

GO 

TO 

400 


?s 


I* ( SC 5 ( «5)) GO TO 500 
N C * fe 

?b IF. C BCS< .6 ) ) GQ T.O .6 0 0 

►i r = 7 

67 IF ( BCSf 7)5 GO TO 600 

NCSA.. . 

?6 !B 60S( 8) ) GO TO 600 

NC = Q 

29 IF ( BCSf. 9)) GQ TO... 600 

NiC = 10 

30 IF ( BCSf 103 ) GO TO 600 

# * IJ s « M P A 

IF ( BC A S ( i ) ) GO TO 700 

32 IP t PC A S f ?) ) GO TO 800 

33 IF f BC AS f 3) ) GQ TO , 900 

C 3 41 

3<J Jp ( BC A S ( <j)) GO TO 1000 

.. .wCss 

35 IP f H C « S \ 5)) GO TO 10 OO 

U C s k 

I F f 6.C A S ( .8) ) . G 0 TO .1 0 0 0 
*C = 7 

JF {BCASf 7)) GO TO 1000 
*C*fi 

TF ( BC A S f A)) GO TO J 0 0 0 
**C *9 

IF fBCASf 0 ) ) GO TO 1000 
NC’sl 0 

IF (BCftS(lO)) GO TO 1000 
IF f ,MOT t STpRFj HETURN 
F** FIL.F <J A 
a&ItF (U 6 f 6013) 

«E TURm 


C.,.,, ST«RT DETERMINING OCMn^aSH CASES, 

r 

■ r 

100 CONTINUE 




c 

C.'./., UNIFORM OOMNMASH 

DO no Hal , M.tfP .... 

DO 110 Psl.FP 
ALPA(P,M)«1 . 

no continue - 

IF (PRINT A ) WRIT E(U6 #60001 
ASSIGN 22 TO NSTaT 

GO TD .2 0.0 0 

200 CONTINUE 
C 

PITCHING ..MODE - - 

!F (CONV?) wRITEfUcO#600l ) 
Rg.4D<U5»5000l X S I C M 

_..DQ .. 2 iC. HeI.tfwp. . 

J*JRATI0*NINDEX(H) 

DO 210 Psl.PP 

xSlejfSlUlPfJ) + CHIC P_(. P J i C Q R D I P. ( J I 
ALFA(P # N)sfXSI/6«AT.XSICN)/C8ARBP 
210 CONTINUE 

I F f PRI nT A ) WRITE fUfc(.6002l XSICm. 
ASSIGN 21 TO nSTaT 
GO TO 2000 

300 CONTINUE ... ... ... . . 

r 

c.,... linear. symmetric twist mode 
00 5t0 P »I# p P 

no 3io h*i.*mp 

AtF A ( P # h ) E ABS ( E T ACP ( M ) ) /8R A T 

310 CONTINUE 

IFfPRINT*) WRITE { l 1 fe # 6 0 0 3 ) 

ASSIGN g4 TO vSTAT 
GO TO 2000 
At 0 0 CONTINUE 

c 

c PAR.A80UC# symmetric T«IST node 

00 <41 0 Msl.MMP 

A|.F a ffTACRfM)/PRAn**2 
00 410 pel, PR 



6I_F6 f P, tf ) = ALF 

aio continue 

IF C PR I NT A ) hRITt(Ufe,600Ul 
ft $ R I n n ?S TO NST ft T 
CO TO 2000 

500 CONTINUE. . - 

C 

r CONTROL SURFACE oomwash mode 

FLPCALs. FALSE., - - — 

.J.F f R R TNT 6, ANC.SYMH5 WRjTE (U6f 6005) 

IF rP»IMTB,A.ND, , NOT, SVMM)WRIT|£fiJ 6,6006) 

ASSIGN ?6 TD. nSTA.T ^ 

GO TO 2010 
60 0 CONTINUE 

r . ,V.'iws^*su>PL"ii!5V nodes 

BJCS(WC)b', FALSE. 

C A L L .. Si QP.E S tUS .. U6# UC0.*.C ONY2, ETA t P *C HlC ?-»-*« — 

lPP»AlFAj- 

IFfPRINTB, aNO.SYmMj KiK T T E ( U6 , 60 0 7 ) 

. .IMPRINT B...AND - .NjD.T..jy.Mi!U .. WRIIE4.U6 j.6 Q.Q.ei_: 

Jr ASSIGN 26 TO M<ST4T 

•° GO TO 2010 

. 700 CONTINUE. 

r 

r, #i , .UN I, FORM aNT-I-SVMMETRIC MODE 

DO .710 Mal-^HMPA 

M-P = 0. 

I F ( E T A C P ( N ) • 6 T , 1 , E ■ 5 1 ALFst, 

. JLEl eta C PJ.N ) , L T . -.1 , E - 51 .ft L F.5 -.1 ^ 

DO 710 P*1,PP 
ALFA f P r M) *ALF 

7J A CONTINUE . . 

IF f PRINT ft) wRiTEfUtf 60095 
ASSIGN 3? TO N S T ft T 

JS Cl , TO 20 00.. _1 . ■ 

600 CONTINUE 
C 

•. . . 



00 PI 0 y-=l , y*PA 
00 MO P=1,PP 
ALFAfP.F isn ACP(m)/BRAT 
Pin CONTINUE: 

TF fP«U'TA) fcPlTEfUM6010) 
aSS’IG^' 3 3 TO ajST AT 
0 TO ? 0 0 0 
non 0 DM T Msfr. 
r 

C *MISYMVFT»I C CONTROL SURFACE DOWNWaSH m 0 Q E ' 

C IFfFLPCfL) CALL FLPDivfur . . , . ) 

I F APRINIP) ^RlTfcrUfef 601 1) 

ASSIGN 3^ TO mSTaT 
gG TO PoiO 
1000 CONTIMiJt 
r. 

C USER SUFPLJE.n NODES 

«CAS(Mc)s. FALSE, 

'■AIL SLOPfcS (Ob# Ofcip 0^0, fOMVP, eta^P, rwirP# mmPa , 

1 PR, ALP A) 

TF (PRIWTB.A^O.SYyy-j WRITE CU6, 60 12) 

TP (PRJ^T H. A«Jn, ,NC)T,SYNM) WRITF(Ufc,fc008) 

° ASSIGN 34 TO Ns T aT 

«0 TO go 10 

r 

C. , , . .POINTING AND STORING MODES 
C 

?000 TFf f NOT ( P*INTi) GO To 2040 
GO TO 2o?0 

2010 IF(.W0T,PRIM9) GO Tn 20410 
2020 r,C 20 3fi «s 1 , mmU 

2030 *PITF(UO,500n ( ALP A t P , M ) , P- j , Pp , 

? 0 4 0 IFfSTORfc.) WWITF.(Ue) f (ALFA(P,M),P ei ,PPi,M- 1# MMU) 
GO TO N5T*T» (22,23,24,25,2b, 32,33, 3^5 
C 

c.,.. .control should never, .get to here 

r 

C,..,. INPUT FORMATS 

5000 FORMAT fgp jo , 0 .) 



500i FORMAT C/riX»9Fl5.7) ) 

r 

. nutp.ut. form a.is~,. 

6000 PORfn&T (V///1BH UNIFORM DOWNfc.ASH'/) 

6001 FORMAT ( 2 2 h ENTER c»*. P0 S,/rREF /> 

6 0.0 2 F Q »r A !..( /.///14H PITCHING KO.Qt . // . 

r 15H C,M, P0S./BHEF.5X, 1H = ,FU.6) 

6003 fORmaT ( / / / / 2 5 h LJNfAR, SYMMfTRIC TNIST /) 

. 60 OR F.P« M A ////£ 5H...P1.R A.BQLI C , _S Y M m.E IRI C .. I W IS T. _V4_._ 'll 

6005 FORMAT ( ////?! H SYMMETRIC f LAP MODE /) 

6006 F0RmAT(////23h UnSYmmETRIC FlAP MODE /) 

_ 60 0.1 FPIRm A T (_/ /Z/.Z.0.H.JJ $. E R»SUPPLIEft SYMMETR I C J3.0.D E / ) . . 

6008 F0PmaT?//// 3?H USER-Si'PPUEO UNSVmmeTRIc mCDE /) 

600 Q F0RMAT(////2OH ANTt-SYMMFTRiC UNIFORM MODE /) 

... 6 0.1 0 — _LCIP.HjA JIjOL L liiS, . JlQbjE _. 

6011 F0PMATf////2fc* aNTJ. SYMMETRIC FLAP MODE / ) 

6012 fOPmaT (///./35H USeR"SUPPL LED ANT I-S YMMETR Ie MODE /) 

6013 r PMA T r////3QR MODES STORED and . F I L E C LOSE D /./ // . 

^UrLOUT 1 N fc. StRAT2 ( XIN,FIN f NIN, XOUT,FOUT # NOUT,NPITE) 

V* C THIS is a' sueKoun ne' to 'determine' the - func t ion Four AT the 

M C XOUT 'LOCATIONS USING UNEaR INTERPOLATION from the 

C ( XI N ,? IN) T A.0L.E , _ : ; _ : . 

r ■ 1 

r XIN MUST EITHER HE IN aScENDING OR DESCENDING ORDER 

. C . XOUT m a .y. be jn,. any nu me rical order’. * 

C N T N- may BE 1, IN whICm Case TnE value F | m [ i ) jS assigned to t f ou T ) a 

c . 

Rfc ALL IN IT 

dimension XIN(NIN) ,EIN(Nl*),*OUmOUT)»FOUT(NiOUT) 

i IN IT ( X > X 1 > X2 , f 1 f pg ) sf 1 ♦ f FP-F 11 / ( X3 *X 1 ) » f X»X 1 ) .... ...... 

C 

c 

IF t NIL.FQ, 1 ) GO TO 600 — 

signsi; .... 

■JFfxiNf i) .GT,xiNf2))siGN-.i T 1 ..' 



NUPRsNJN.I 

DO 500 NO*l f NQUT 

l£AS.l S.N»XOUTrNQ) t r,T.SIr;N»XINf21 ) r.O TP 10Q 

FOUTCNO) c LXNIT(X0UTCN0)#XIN(i)»XlN(2) f piN{i)#rlN<2)) 

GO TO 500 

— 1-0-0 IPjI1SIG m»XQUT f nQ) .LT.S lG NfcXlNf uUPRn r.n ip 2on _ 

FCUTfNO) = LINJTfXOUT(NO) ,XIN(NUPR) , XJN (NINJ , FIN (NUPR) ,FIN(NIN) } 
GO TO 500 

POO CQ.NT IN UE: 

DO 300 NI*3| NUPR 

lF(SIGN*XDUT(NO) , G? , S I GN*X IN f N I) ) GO TO 300 

F OU T t N O. ) a LTNTTfXOUTfNOi.XlNtMi-n f X;N(Nt1 f 

lFlN(NI-n #FlN£NI) ) 

GO TO 500 

300 CONTINUE 

C 

WRITE (NR I Tp, 1 ) XXN 

— 1 FORM 4 T' ( t Hj t L IABLE NOT IN ASCENDING OR DESCENDING ORnE R I Mi 

l f ' SUB, 5TRAT2'/(lXf lPp2e t 5)) 

STOP 13 

C 

wi 500 CONTINUE 

M return 

(l£LG DO 610 NOcI.nOUT 

610 fOUT(NO)*FIN'( 1) 

RETURN 

END ' 

RUgROUTlNE SLOPES (U5,U6,UC0,CGNV2 ,ETaCP. CHICP.MMU, 
lPP#ALFA) 

integer us.u6.ucn,pp f p _ 

LOGICAL CQNV2,L0G1,TE$T 
REAL L 

DIMENSION ETaCR(NMIj) r CHICPrPP), ALF A fPP r MMLH 

DIMENSION eTA(20),ALfA 2(15,20), ALFA3(30), cHlfai)* 
lALF(4l)f Af40), B(flO), C(«0), L(40), G(<I0), 

miO-If H f AfO 3 

DIMENSION TEST (15)#CHIDUM(3Q) 

C 

■_C^ ^.THE^IME_NSIQN OF THE ABOVE TWO arrays SHOULD BE AS LARGE 



nnpnnh-nnh pm non I o o p n hnnbnn 


c , , . . . 

c 

C » n i i 

c • « • f • 



AS PP AND 2*PP . RESPECTIVELEY 

J HE DIMENSION OF ALFa^ MUST flF THF MaXTMHM nF p „PP A Np 

nsmax, 

NSmaX is ThE MAXIMUM NUMBER OF SPaNWiSE STATIONS at 

whi c h d a t a sho uld r e giv en '. 

NSMAX is A DIMENSION FOR ETa,ALFA2, and alfas« 

Data NSMAX /20/ 


...'.•ntbmaX.i is the maximum num 9 er of chordwise locations 
.....at which the camber or slope data should be given, 

...... ntbma x is T ^L_DIM.£_NSi Q_N_.QF cm and alf* . 

, , , , .ntbmax-i is the dimension of a, b, c, l, g, 
m.mEi aND H. SINCE a is also used for a camber or SLOPE 

.....POLYNOMIAL. ThE DIMENSION OF a MUST BE AT (_FAST B, 
DATA NTBmax /«■! / 


I_ PS I S USE D F OR numerical differentiation 

OF THE CaMBER DISTRIBUTION DETERMINED BY A TABLE AND 
COOI^g (CONTROLLED DEVIATION INTERPOLATION METHOD), 
D ATA EpSPLF /S .E f. 3-/ 


NPPMaX is THE FIRST DIMENSION. Of aLFA2 and THE MAXIMUM 

ALLOWABLE VALUE FOR PP‘ t _ 

DATA NPPmAX /15/ 


IF(PP # GT # NPPMAX) Call ST0P2(U6,» PP TOO LARGE IN SUB, SLOPES I, 
IFLOaT(PP)) 


...,.NSPSEC is THE NUMBER OF SPaNWISE SECTIONS 
.....ITYPES DENOTES THE TYPE OF SP aN WISE INT ERPOLATION, 
..... 0 IMPLIES STRAIGHT line interpolation’, 

..... NOT 0 IMPLIES THAT C00IM2 WILL BE USED' 

IFfcQNV?) WRlTEfUfOffcnftnl 


ij *j u u u u d t-j «J u u uu u ci c: m «J» ts 

> < . i ! • i 


I F f NSP$£C".GT.N$MAX) CALL ST 0 P 2 fU< 5 # 

1» TOO MANY SPAN Wist SECTIONS In SUB , SlOPES, • * FLOAT CnSPSEC) ) 

RE A D.tU.5 # 55. 0 G.)_N 3 P SE Cj I.T YP.ES • 

IPf ITY°tS,FO.O) GO TO j>0 
JF CC0NV2) WRITE rUCO. 6001) 

9EA0.(U5,5fl.01i-XK. 

?0 CONTINUE 
c 

QQ._bL.Q _0 SPS EC. 

IF'fCONV 2 )WRITE(Uc 0.6 006 ) 

REAOft'S#5001) FTA(N), TWIST 

c . • 4“ f$ s p7 n w‘i sFTa r ' i a ble~ it ~m7V 'exceed i'. in value 

C • .... IN ORDER TO CONTROL THE c00IH 2 INTERPOLATION NE aR THE ENDS, 

W1SI-13. J HE. A. NGLE--iiF ..TWIST,. - • • - 

IPfC 0 NV 2 ) WRITE (UCO# 6002 ) 

RE AO (U5# 5000 ) ITYPEC 


ITYPEC DENOTES The TYPE OF c^ORdWISE SLOPE DEFINITION 
THAT the USER wants to use 1 , 

0 TMPL IES .That JTiiERE jS.m.C> M .BER,.„jUS.T . TWIST. .... - 

1 IMPLIES That A set of polynomials DEFINES the slope 


DISTRIBUTION 

,, , , ^ 2 implie s Tha t a s et of polynomial S..OE f I N E S_.ThE 

Camber DISTRIBUTION', 

3 implies that a table of values defines the slope 

; ; ; ; ; DI STRIrU TION SU rRQUT I_NE COD JMg _w.IL L_ S E USEp f OR .THE. 

. . .7. interpolation^ 

..... u ihplifs that a table of values defines the 

1JL / tA Camber D is T RIrUTICn: c OP Im^. WILL _ aE . y SfJL T o 

. determine tmE camber in The vicinity of The control 
..... POINTS, THEN NUMERICAL DIFFERENTIATION WILL RE USED 

.. ... TO DE TE RMINE thp SLQPFS. 

,,,,, 5 IMPLIES TnAT A TABLE OF VALUES DEFIES ThE SLOPE 

DISTRIBUTION, A CUBIC SPLINE FIT WILL BE USED TO 

..... oftfr mine values at the c ho r p his e _g ont _rol.. st aIion s . 

..... 6 IMPLIES That A TARlE OF VALUES DEFINES ThE CamBEr 

..... DISTRIBUTION, a CUBIC SPLINE WILL BE DETERMINED FOR 

..... THE Tab L E and then DIFFERE NT I A IEJEJLQJO rTjI I.N. the 



c 

e 


SLOPES, 


— If_rlTYPEr.EQ.oi aft Th ion 

IPCITYPEr.LT.O) GO TO a50 
lFflTYPEc.LT. 31 GO TO 200 

ULai-YPEC.LT.T) GO TO 300 

GO TO ii 5 0 
100 CONTINUE 

D..Q 110 P*1,PP 

110 ALFA?(P,N)s TWIST 

GO TO 500 


2.0.5 CQNT JLNUE. 

C 

C ' 


C,., . .EVALUATING CHORDWISE SLOPES FROM polynomials; 

c 

- IE^CQNV2) -WRIXlf Uca# 6003} 

PEADfUs^sooO) NPOl S 
CwIMINsO. 


U1 

l-n 


cm.. 210 p=nei? 

210 TEST (P)x, TPUfc', 

no 280 NP5_1*NP0LS 


60.0J4_ 


C 

C CHIMAK DENOTES THE UPPER LIMIT OF THE POLYNOMIAL THE 

c***.**. default value is i, * 

C SCAiF IS A SCALE FACTOR FpR THE POLYNOMIAL 

C COEFFICIENTS, DEFAULT Y A LUE is \' m 

C» . ■ « . C Hi ref is the ORI GIN OF THE polynomial; _ 

9.F4DruS,500n ChIMAX, SCALE»CmI 9EF 
IFfCHlHAX.EQ.O.J CHlMAXsl, 

RF.APfU5> 500l ? (Aflinsiffll 

c 

C THE number of COEFFICIENTS IS LIMITED TO e; 

C ...... , , THE F OLLOW ING DETER MINES T HE ACTUAL N UMRER’ 

NT ERwSs0 

DO 2 1 5 1st, 8 

IF.XAI_NTXRM51^NE^.5*.1 BO_.LO._225 



2 l 5 


220 

225 

226 
C. 

C * • « .« 
(!•••< 


230 

240 


?45 

270 

275 

?eo 


285 


NTFRmSsNTERmS-1 

CONTINUE 

CALL 5T0P2(U6, •. .A.CAHBER. OR SLOPE POLYNOMIAL. IN SUB, SLOPES 
1 has ALL 0 COEFFICIENT fTAs t*f.TA(N)) 

IFfSCALE.EQ.O. )GO TO 22#, 

DQ 225 I=1,NTERmS 

A(n*AcI)*8CALE 

CONTINUE 

the following evaluates the slope for all control points 
’which LIE RIThIN the Range of thf polynomial. 

DO 275 psi.,pp 

IF CCHlCPfP) ,LT,CHI m Jn) GO TO 275 
IF ( C M I C’P t P ) .GT.CMI^aV) GO TO ? 75 
IF f ITYPfcc.EU.2l GO TO ?40 

A | FA?(P,N)rAf t ) - T * T S T 

I f (NTFRMS.EG.n GO TO 27 0 
YsCHlCP(P)*CHlRtF . 

DO 23n 1x2, N TERMS 

ALrA?(P,N)sALFA?fP,N) *• X*A( j) 

XsX*(CMICP(P)-CMIREF 1 
GO TO ^ 7 0 

A| FA2(PfM)=A(2) - T'MST 
IF ( NTfRMS , L T , i ) GO TO 270 

XsCHlCPfPUCHlREF 
00 245 I s 3 , N T f R M S 

ALFA?(P»N)sAL FA gfP.N1 + FLOAT (T„i)*X* 4(1) 
XxX*CCHlCPtP)-CH|RFF > 
continue 
tfst(P)=, false, 

CONTINUE 
CHI KTN-CMlMAy 

continue 

LOGU, FALSE. 

DP 285 Pxi.PP 
LOGlsLOGl .OR.TFST (») 

IFfLOGl) till ST0PPCU6- , _ , 

p a CHhrDwIsE PnLY, OIO n ,) T COVER ENTIRE CHnRD. r.T A s «, 

2ETa ( W)) 



GO TO 500 

r 

t evaluating the chorpwise slopes from” Tables 

300 CONTINUE 

I f it 0 N V 2 ) _W R I T E C U CO i 6 . 0.0 5 3 

Rp;AD fU5, 500 1 ) $Calf>A 1,AN 
¥«CeAl 

C SCALE IS A SCALE f 4 L TOR FOP THE UBLE ENTRIES 

* C,» , ( , a 1 AND AN a«E THE SLOPFS A T THE END PQImTS Of'thF T 4 RLF 

■C « A 1 _A N. D . . ..AN ,._S HO OLD. BE P U T .. I N . 1 F PQS.SIRLE j ESPECIALLY IF * 

c , . , , , T he table IS one giving the Ca m BEH. IF both VALUES are 

C ACTUALLY ZERO,. THEN INPUT ONE AS 1,6-30, IF EITHER IS 

— fU* JSIY£±U_J tithi ._.&£) J -H MUST. b£_JS1VD4 _. . 

C. f ,.,**e IS THE end POINT INTERPOLATION CONTROL FOR c00lM2 
C.....IF XKC*0«, Then STRAIGHT LINES will be used IN the end 

— — j.£_KKC.«.Lm: _ IH fLN F U L l PARABOLIC 

c. .... INTERPOLATION hill be USED in The end INTERVALS a 

c value in between hill give a curve in between. 

: -lliS-CALS.. £ Q . 0 . ) SCALEfcl - . .... 

\n *\ f A 8 s 0 

-' J 310 CONTINUE 

N T A B g N T A B « i 

REAorUBf 5001 ) CHI(NTA 8 )f ALF(NTAB) 

IFfNTABj.GT.NTBHAXi CALL STQP? f Ufe, 

™L...... 1 L::.IOO..NaNY TA.BLE._FN TRIES I N SUg. SLOPES. £T A* «,ETA(N>) 

IF(CmI(NTAB> .LT.98, ) GO TO 310 
NTAB 3 NTA 8 -I 

I F..CI T.YPE.C . EO , 3.3 GO_ T 0 315 

lF(I-TYPEc«ea.fl) GO TO 350 
GO TO «00 
315 CONTINUE 

c 

c the table is one of slopes and codim 2 will be used 

GALL ..cop Mg JC aL f f c H 1 1 NT A &J..C HICPi ALFA 3 * PP» XKC ) 

DO 320 Psl,PP 
ALFAg(P # N)sALFA3rP)*THIST 
_ S2fl cMIIttUE ; ...: 



r,C TO 500 
35 ft CONTINUE 

C . . . . 

C THE TaBLE IS ONE OF CiMjJEW and C00IH2 a ND numerical DIF* 

r feremttat ion ^ t l l used, 

Js-l 

C 

r F.o SDEf IS THE DEFAULT VALUE FDR EPS (SEE DATA STATEMENTS) 

C • ... - : 

FPS s EPSOEF 

IF ( AM ,nE , 0 , JEPSsAM 

DC 355 Psl.PP 

J« J+2 

CWlDiJMf j)sCHiCP(P)-EPS 
CHlDU«CJU)sCHICPCP)+EPS 
355- COMTINUF 

CaLL COOIMsc aLF,chI,nTa8 ,CHI0U'M, ALFA3#a*PP,XKC) 

•Js*l . . 

On 360 Fsl.PP 
J*J + 2 

360 ALFA2(P,N)«tALFA3{Jtn-ALFi3(J3)/(2.*EPS) - TRIST 

GO T n 500 
UOO CONTINUE 

C . 

CUBIC SPLINE FIT WILL bE USED TO DETERMINE THE CaMBER 

OR SLOPE DISTRIBUTION, SPLJNj IS CALLED 

ThE End point DERIVATIVES are NOT GIVEN. IT sets 

C The SECOND DERIVATIVE to zero at the end POINTS, if THf 

c E^'D point DERIVATIVES are KNOWN to a high decree of 

c. .... accuracy, c,i .percent?) TmEn SPL I N2 ShOUlD BE USED AS jT 

C 'V 0 R K S MUCH BETTER. 

IFfAl ,FG,o, , AND, AN,EO,fi,l 

1 C A LL . INI (NT aB|CuI, AL.EfA f B f CfL#D»^^H) 

IFfAl, NE,0, , OR , an,nE,0,i 

ICaLL SPLIN2(NTaB,CHI , *LF, A , B , C , L , g , E , H , A 1 , ANl 
c 

C IN The above a,B,C are The local SPLINE COEFFICIENTS 

C AND L»G#F, and H are WORKING SPaC'V. 

do «30 Psl,PP 



DC 410 Jsl , NTAB 
MsJ.I 

QlX* C-Hl - CfeLLf. J.) 

IPfDX.LT.O.) GO TO 40 
410 CONTINUE 

MlsNTA.9 . .. 

415 Nl-«MAXi>(l#Nl) 

D*=ChICP(P (Ml) 

- I F 1 1 T V p L C.* E@ . 6 ). . . G. 0 .. JjQ „.42 a 

ALPA?(P#N)-ALE(Nl)4'OX*(Afrvl)+OX*fBfNl)'*DX*CCNl)))-TWIST 
00 T n a in 

QM 4 i E A ? ( P j N ) s A f M.1 1*2 , N 11 *&X * 3 f » C ( N 1 ) * 0 X * * 2 - TWIST • 

C 

C.,... IN EVALUATING the above SPLINE CURVES TT Has ASSUMED that 

f. •■*• ».» N O CONTROL’ P.0 1 .NX5 _LXE.D..*.QlUJ.3.1D£ 0 E. ..T.H£._J.a B1 E, . . 

«3rt CONTINUE 

50 TO 500 

455 CONTINUE ; ; 

C 

c . ... .evaluating the slopes using a user-supplied 
c . . ..... sub.r.out .1 Mt, :r.„. .. . ; • 

ui r 

J Call ■US-U>PE(U5,U6fUC0»C0NV2 # N.ETAfN),PP,CHlcP, 

1 .1 T Y ° E C # a L R A 3 ) 

bo -460 Psl ,RP 
ALFA?(P,N)-ALKA3(P1+TWIST 

4"6o_.... CONTINUE... 

500 CONTINUE 

r 

r 

t .. . - - 

C 

c ii... . : _ 

c ' ■ " '■ 

C.,...a t This POINT tm F DETERMINATION Or. the SLOPE DISTRIBUTION 
C.....£ t the INTERSECTIONS OF the CHORDWISE CONTROL LINES AND 
r..,,, T MP GIVEN eTa STATIONS Has BEEN ACCOMPLISHED. NOW it 
C.....TS NECESSARY TO INTERPOLATE TO OBTAIN ThE SLOPES 
C. . . . • 1 T. ,.T H Ej S P A N w i s £ C 0 N T R 0 L . S T a T I Q n s . G I Vj* N . BY. . E T A C.° . 



r 


00 1000 P=1,P» 

DO .600 .N..S.1 NSP.SE C 

c 

C * • THE minus sign BELOW OCCURS BECaUSE THE INDUCED DOWnWaSH 
C.... .EQUALS minus THE given SLOPE*. 

*L FA ''(N)s-4 L FA2fP,N) 

600 CONTINUE 

I.F £lHPES,EQ. o). .CULL STRa T 2 £ E T A # A LF.A l # N SPSEC # E T aCR # 

1 ALF# HMU,UCO) 

IFfITYPES.NE.Oi CALL CODIMpf aLFa3, ETa .NSPSEC # STaCP, 

. lALF# MhU#XK3 

DO 6?0 

AL F A(P.M)sALPfM) 

_62fl C.CN.T INUE 

1000 CONTINUE “ 

Rg TUR V 

-C- 

C 

C.... .formats 

c .. . 

2 5000 FORMAT (1615" ) 

50 01 FORMAT (BF 10,0 ) 

6000 FpRM4T(i ENTER .NSPSEr.Si ITYPES*/) 

6001 FOPm a T f i ENTER X K i / ) 

600 ? FORM 4 T ( » ENTER I TYPED 1 /) 

6003 PQRv AT ( « EnTE.R nPOl s ' f ) 

SO 0 a E 0 R M A T ( i ENTER CHI* aX , SC ALE , CHI»EF> t / i POLYNOMI AL ’cOEFFICIENTST/) 

600s FORMAT C» ENTER TaBLE SCALE FACTOR# aICOR XKC)#f/ 

1' A^COR EPS)# And THE TARLe. iC.HI=99. STOPS) »/) 

6006 EORm a T f i ENTER ETa, T*I5Ti/) 

Subroutine $p_l I N l ( n , * 7F # A , R # C # L , G ~ £_, h j 
REAL L 

dimension v(n) 

DIMENSION F(.v) 
nlMFNSIUN A (N) 

DIMENSION b f w j 
DIMENSION £(W) 


DIMENSION L(N) 

DIMENSION G(N) 
n I fc N s I C N _ e Cisi-3- - - 

n I ME^SlON H(N) 
y. r N m 1 

C ( -M ) s 0 , „ __ 

Bfl)sO, 

B.( w .) a 0 « 

- - 00.10 ,1 si r M 

i.r !)sxc i + n-xn) 

G(!) = {F(I + 1 )-F(n )/L ( 1) 

1.0 CONTINUE. ... . 

E ( 1 ) sO , 

DO 20 1*2, « 

E(I)»Lfh/D 

- H ( I ) s.( X I.) - G ( J ) 1 -L ( J).) /.D 

2o CONTINUE 
FfMj-C'. 

2 00 30 ,J*2,K 

I)*8f 1 + 1) 

IfleJ 

30 CONTINUE ; 
do ao isj,M 

JLUJs. G.UJ CLL1)_* a8 £1 It B ( I tl )..) / 3 , 

C(I)s{BcI*n-B{n)/C3,*LtD) 

4n CONTINUE 

_ i LNJ.tA_(_« ) + L_£ H. I* B. f M ). 

Rf TURN 

END ■ 

BUflROUTi.N| . 5?.LLH2xN*XJf * Ju.B ,.C,.L , AN ) 

NEaL L 

D I NEWS I ON X(N) 

MmENSIO.n F(W) 

DIMENSION A ( N J 
DIMENSION g(Nj 

DIME NSION rfM . 



0 I HENS I Ok L(M 
MANSION G ( N) 

Ql PENSION ECNl 
M ME m$ I ON H C Ki > 
s N - 1 

00 10 Isl, M 

G(l)s(F(I+1 )«F(I))/L(I) 

10 CONTINUE 

Del ,5*1, ( 1 ) +2 ,*L f 2 ) 

EC2)sLr2)/0 

H (2)*.l3.*G{2) - a,5*G C n + i.b*All/0 

DO ?r- I = 3 P M 
Jsl-l 

Os? . * ( L c J J *1 ( n ) -l { J ) * E ( • J ) 

E’fI)«Uf I)/9 

u fI)*C5.*(G(I)-G(«J) J)}/D 

jr rO^TTNUE 

OsJ,*Lff1) + 1 ,5*1. f M+i )-Lfri)*FfM) 

p fM+!. ) = 1,5*4 

** DC 3 0 J = 

TsT.j 

3n r 0 Ni T I MJ E 

D(l)=3',*(G(i)-4l)/f?;*ua) > - ',5*HC2) 

msM-1 

Ff 1 )r3,* (AM-G(MJ 
00 UC I = 1 f M 

A(l)sG(!l»U (!)*(?. *e(I)+w f l + i))/3. 

CfI)*(B(! + n - Bfm/( 1 .*L(!)) 

Ur CONTINUE 
4 ( S! ) s 4 w 
C(E')sO, 

Rf,TU«N 



SUBROUTINE ST0P2fN,HES46E,V A t ) 

DIMENSION v t S4GE r?r ) 


WHITE (N»J) ME5*6EfV*L 

STOP n 

1 FDh’y.&I (///*. ***** i,2Q.A(t,t *****!/! VA(.,s I , i PE IS , 7 ) 
E f r> 



